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We consider quantum Hamiltonians of the form H = Hq — U^. cos{Cj) where Hq is a quadratic 
function of position and momentum variables {xi,pi, X2,P2, ■■■} and the Cj’s are linear in these 
variables. We allow Hq and Cj to be completely general with only two restrictions: we require that 
(1) the CjA are linearly independent and (2) \Cj,Ck\ is an integer multiple of 27ri for all j,k so 
that the different cosine terms commute with one another. Our main result is a recipe for solving 
these Hamiltonians and obtaining their exact low energy spectrum in the limit 17 —>■ oo. This recipe 
involves constructing creation and annihilation operators and is similar in spirit to the procedure 
for diagonalizing quadratic Hamiltonians. In addition to our exact solution in the infinite U limit, 
we also discuss how to analyze these systems when U is large but finite. Our results are relevant to 
a number of different physical systems, but one of the most natural applications is to understanding 
the effects of electron scattering on quantum Hall edge modes. To demonstrate this application, we 
use our formalism to solve a toy model for a fractional quantum spin Hall edge with different types 
of impurities. 


I. INTRODUCTION 

In this paper we study a general class of quantum 
Hamiltonians which are relevant to a number of dif¬ 
ferent physical systems. The Hamiltonians we con¬ 
sider are defined on a 2N dimensional phase space 
{xi, ...,xn,Pi, ■■■,Pn} with [xi,pj] =idij. They take the 
form 

M 

H = Ho-UY,cos{a) ( 1 ) 

i=l 

where Hq is a quadratic function of the phase space vari¬ 
ables {xi, ...,Xn,Pi, ■■■,Pn}, and Ci is linear in these 
variables. The Cfs can be chosen arbitrarily except for 
two restrictions: 

1. {Cl, ...,Cm} are linearly independent. 

2. \Ci, Cj] is an integer multiple of 27ri for all i,j. 

Here, the significance of the second condition is that it 
guarantees that the cosine terms commute with one an¬ 
other: [cos(C'j),cos(C'j)] = 0 for all i,j. 

For small U, we can straightforwardly analyze these 
Hamiltonians by treating the cosine terms as perturba¬ 
tions to Hq. But how can we study these systems when U 
is large? The most obvious approach is to expand around 
U = oo, just as in the small U case we expand around 
[7 = 0. But in order to make such an expansion, we first 
need to be able to solve these Hamiltonians exactly in the 
infinite U limit. The purpose of this paper is to describe 
a systematic procedure for obtaining such a solution, at 
least at low energies. 

The basic idea underlying our solution is that when [7 
is very large, the cosine terms act as constraints at low 
energies. Thus, the low energy spectrum of H can be de¬ 
scribed by an effective Hamiltonian i7eff defined within 


a constrained Hilbert space TLeS- This effective Hamilto¬ 
nian Heff is quadratic in {xi, ...,xn,P i, ■■■,Pn{ since Hq 
is quadratic and the constraints are linear in these vari¬ 
ables. We can therefore diagonalize Heff and in this way 
determine the low energy properties of H. 

Our main result is a general recipe for hnding the ex¬ 
act low energy spectrum of H in the limit [7 —>• oo. This 
recipe consists of two steps and is only slightly more com¬ 
plicated than what is required to solve a conventional 
quadratic Hamiltonian. The first step involves finding 
creation and annihilation operators for the low energy 
effective Hamiltonian i7eff (Eq. 11). The second step 
of the recipe involves finding integer linear combinations 
of the Cfs that have simple commutation relations with 
one another. In practice, this step amounts to finding 
a change of basis that puts a particular integer skew- 
symmetric matrix into canonical form (Eq. 14). Once 
these two steps are completed, the low energy spectrum 
can be written down immediately (Eq. 19). 

In addition to our exact solution in the infinite [7 limit, 
we also discuss how to analyze these systems when [7 is 
large but finite. In particular, we show that in the finite [7 
case, we need to add small (non-quadratic) corrections to 
the effective Hamiltonian Heff in order to reproduce the 
low energy physics of H. One of our key results is a dis¬ 
cussion of the general form of these hnite U corrections, 
and how they scale with U. 

Our results are useful because there are a number of 
physical systems where one needs to understand the effect 
of cosine-like interactions on a quadratic Hamiltonian. 
An important class of examples are the edges of Abelian 
fractional quantum Hall (FQH) states. Previously it has 
been argued that a general Abelian FQH edge can be 
modeled as collection of p chiral Luttinger liquids with 
Hamiltonian^^"^ 

i7o = J dxj^{d,<4>fV{d,<I>) 


2 


Here $(a;) = {(j)i{x), ...,(j)p{x)), with each component 
4>i describing a different (bosonized) edge mode while 
V is a p X p matrix that describes the velocities and 
density-density interactions between the different edge 
modes. The commutation relations for the (/>/ operators 
are [(j)i{x),dy(l)j{y)] = 2Tri{K~^)jjS{x — y) where if is a 
symmetric, integer p x p matrix which is determined by 
the bulk FQH state. 

The above Hamiltonian Hq is quadratic and hence ex¬ 
actly soluble, but in many cases it is unrealistic because 
it describes an edge in which electrons do not scatter be¬ 
tween the different edge modes. In order to incorporate 
scattering into the model, we need to add terms to the 
Hamiltonian of the form 

Hscat = j U[x) cos{K^— a{x))dx 

where A is a p-component integer vector that describes 
the number of electrons scattered from each edge mode®. 
However, it is usually difficult to analyze the effect of 
these cosine terms beyond the small U limit where per¬ 
turbative techniques are applicable. (An important ex¬ 
ception is when A is a null vector®, i.e. KK = 0: in 
this case, the fate of the edge modes can be determined 
by mapping the system onto a Sine-Gordon modeh). 

Our results can provide insight into this class of sys¬ 
tems because they allow us to construct exactly sol¬ 
uble toy models that capture the effect of electron 
scattering at the edge. Such toy models can be ob¬ 
tained by replacing the above continuum scattering term 
Hsc&t by a collection of discrete impurity scatterers, 
cos(A’^if$(xi) — Ui), and then taking the limit 

—>■ oo. It is not hard to see that the latter cosine terms 
obey conditions (1) and (2) from above, so we can solve 
the resulting models exactly using our general recipe. Im¬ 
portantly, this construction is valid for any choice of A, 
whether or not A is a null vector. 

Although the application to FQH edge states is one 
of the most interesting aspects of our results, our focus 
in this paper is on the general formalism rather than 
the implications for specific physical systems. Therefore, 
we will only present a few simple examples involving a 
fractional quantum spin Hall edge with different types 
of impurities. The primary purpose of these examples is 
to demonstrate how our formalism works rather than to 
obtain novel results. 

We now discuss the relationship with previous work. 
One paper that explores some related ideas is Ref. 8. In 
that paper, Gottesman, Kitaev, and Preskill discussed 
Hamiltonians similar to (1) for the case where the Ci 
operators do not commute, i.e. [Ci,Cj\ ^ 0. They 
showed that these Hamiltonians can have degenerate 
ground states and proposed using these degenerate states 
to realize qubits in continuous variable quantum systems. 

Another line of research that has connections to the 
present work involves the problem of understanding con¬ 
straints in quantum mechanics. In particular, a number 
of previous works have studied the problem of a quan¬ 


tum particle that is constrained to move on a surface by 
a strong confining potential^bo. This problem is similar 
in spirit to one we study here, particularly for the special 
case where [Ci, Cj] = 0: in that case, if we identify Ci as 
position coordinates Xi, then the Hamiltonian (1) can be 
thought of as describing a particle that is constrained to 
move on a periodic array of hyperplanes. 

Our proposal to apply our formalism to FQH edge 
states also has connections to the previous literature. In 
particular, it has long been known that the problem of 
an impurity in a non-chiral Luttinger liquid has a simple 
exact solution in the limit of infinitely strong backscat- 
tering^^^^^. The infinite backscattering limit for a sin¬ 
gle impurity has also been studied for more complicated 
Luttinger liquid systems^®^^®. The advantage of our ap¬ 
proach to these systems is that our methods allow us to 
study not just single impurities but also multiple coher¬ 
ently coupled impurities, and to obtain the full quantum 
dynamics not just transport properties. 

The paper is organized as follows. In section H we sum¬ 
marize our formalism and main results. In section HI we 
illustrate our formalism with some examples involving 
fractional quantum spin Hall edges with impurities. We 
discuss directions for future work in the conclusion. The 
appendices contain the general derivation of our formal¬ 
ism as well as other technical results. 


II. SUMMARY OF RESULTS 

A. Low energy effective theory 

Our first result is that we derive an effective theory 
that describes the low energy spectrum of 

M 

H = Ho - U'^cos{C,) 

i=l 

in the limit U ^ oo. This effective theory consists of an 
effective Hamiltonian Hgg and an effective Hilbert space 
HeS- Conveniently, we find a simple algebraic expression 
for Hgff and HeS that holds in the most general case. 
Specifically, the effective Hamiltonian is given by 

HeS = Ho-J2 2 '' ■ 

where the operators Hi, ...jHm are defined by 
1 

= (3) 

and where Mij is an M x M matrix defined by 

M = [Cj,Ho]] (4) 

This effective Hamiltonian is defined on an effective 
Hilbert space Heff, which is a subspace of the original 
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Hilbert space "H, and consists of all states \il}) that sat¬ 
isfy 

cos(C'i)IV') = |'0), i = (5) 

A few remarks about these formulas: first, notice that 
M. and N are matrices of c-numbers since Hq is quadratic 
and the C^’s are linear combinations of Xj’s and Pj's. 
Also notice that the H^ operators are linear functions 
of {xi, These observations imply that 

the effective Hamiltonian is always quadratic. An¬ 
other important point is that the H^ operators are con¬ 
jugate to the Ci’s: 

= ( 6 ) 

This means that we can think of the H^’s as generalized 
momentum operators. Finally, notice that 

[C„iJeff]=0 (7) 

The significance of this last equation is that it shows that 
the Hamiltonian Hes can be naturally defined within the 
above Hilbert space (5). 

We can motivate this effective theory as follows. First, 
it is natural to expect that the lowest energy states in 
the limit [/ —> oo are those that minimize the cosine 
terms. This leads to the effective Hilbert space given 
in Eq. 5. Second, it is natural to expect that the dy¬ 
namics in the Ci,...,Cm directions freezes out at low 
energies. Hence, the terms that generate this dynamics, 

namely —^H^Hj, should be removed from the ef¬ 

fective Hamiltonian. This leads to Eq. 2. Of course this 
line of reasoning is just an intuitive picture; for a formal 
derivation of the effective theory, we refer the reader to 
appendix A. 

At what energy scale is the above effective theory 
valid? We show that iJeff correctly reprod uces t he en¬ 
ergy spectrum of H for energies less than yjljjm where 
m is the maximum eigenvalue of M.ij. One implication 
of this result is that our effective theory is only valid 
if N is non-degenerate: if Af were degenerate than AA 
would have an infinitely large eigenvalue, which would 
mean that there would be no energy scale below which 
our theory is valid. Physically, the reason that our effec¬ 
tive theory breaks down when Af is degenerate is that in 
this case, the dynamics in the Ci, ...jCm directions does 
not completely freeze out at low energies. 

To see an example of these results, consider a one di¬ 
mensional harmonic oscillator with a cosine term: 

rA Kx‘^ 

H ^ ^ — -U cos(27r:r) ( 8 ) 

2 m 2 

2 If' 2 

In this case, we have Hq = + and C = 27ra:. If we 

substitute these expressions into Eq. 2, a little algebra 
gives 


As for the effective Hilbert space, Eq. 5 tells us that 'Heff 
consists of position eigenstates 

■Heff = {|a; = 9 ), 9 = (integer)} 

If we now diagonalize the effective Hamiltonian within 
the effective Hilbert space, we obtain eigenstates \x = 9 ) 
with energies E = Our basic claim is that these 

eigenstates and energies should match the low energy 
spectrum of H in the 1/ —>■ 00 limit. In appendix Ala, 
we analyze this example in detail and we confirm that 
this claim is correct (up to a constant shift in the energy 
spectrum). 

To see another illustrative example, consider a one di¬ 
mensional harmonic oscillator with two cosine terms, 

Kx"^ 

H = - - 1 --- U cos(dp) — U cos(27ra:) (9) 

2m 2 

where d is a positive integer. This example is fundamen¬ 
tally different from the previous one because the argu¬ 
ments of the cosine do not commute: {x,p\ ^ 0. This 
property leads to some new features, such as degenera¬ 
cies in the low energy spectrum. To find the effective 
theory in this case, we note that Hq — ^ and 

Cl = dp, C 2 = 2ttx. With a little algebra, Eq. 2 gives 

-ffeff = 0 

As for the effective Hilbert space, Eq. 5 tells us that HeS 
consists of all states IV') satisfying 

cos(27rx)|'0) = cos(dp)|'0) = |'0)- 

One can check that there are d linearly independent 
states obeying the above conditions; hence if we diagonal¬ 
ize the effective Hamiltonian within the effective Hilbert 
space, we obtain d exactly degenerate eigenstates with 
energy E = 0. The prediction of our formalism is there¬ 
fore that H has a d-fold ground state degeneracy in the 
17 —>■ 00 limit. In appendix Alb, we analyze this exam¬ 
ple and confirm this prediction. 

B. Diagonalizing the effective theory 

We now move on to discuss our second result, which is 
a recipe for diagonalizing the effective Hamiltonian TIeff- 
Note that this diagonalization procedure is unnecessary 
for the two examples discussed above, since ideff is very 
simple in these cases. However, in general, ideff is a com¬ 
plicated quadratic Hamiltonian which is defined within a 
complicated Hilbert space dfeff, so diagonalization is an 
important issue. In fact, in practice, the results in this 
section are more useful than those in the previous section 
because we will see that we can diagonalize ddgff without 
explicitly evaluating the expression in Eq. 2. 

Our recipe for diagonalizing ddeff has three steps. The 
first step is to find creation and annihilation operators 
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for -ffeff- Formally, this amounts to finding all operators 
a that are linear combinations of {xi, ...,xn,Pi, ■■■,Pn}, 
and satisfy 

[a, = Ea, 

[a,Q]=0, i = (10) 

for some scalar E ^ 0. While the first condition is the 
usual definition of creation and annihilation operators, 
the second condition is less standard; the motivation for 
this condition is that iJgff commutes with Ci (see Eq. 7). 
As a result, we can impose the requirement [a, Ci] = 0 
and we will still have enough quantum numbers to diag¬ 
onalize iJeff since we can use the C^’s in addition to the 
a’s. 

Alternatively, there is another way to find creation and 
annihilation operators which is often more convenient: 
instead of looking for solutions to (10), one can look for 
solutions to 

M 

[a, Hq] = Ea + ^ \j [Cj , H ^], 

[a,Ci]=0, i = l,...,M (11) 

for some scalars E, Xj with E ^ 0. Indeed, we show in 
appendix B 4 that every solution to (10) is also a solution 
to (11) and vice versa, so these two sets of equations are 
equivalent. In practice, it is easier to work with Eq. II 
than Eq. 10 because Eq. 11 is written in terms of Hq, 
and thus it does not require us to work out the expression 
for Heff. 

The solutions to (10), or equivalently (11), can be 
divided into two classes: “annihilation operators” with 
E > 0, and “creation operators” with E < 0. Let 
ai,...,aic denote a complete set of linearly independent 
annihilation operators. We will denote the correspond¬ 
ing E’s by Ei,...,Ek and the creation operators by 
a|,...,a^. The creation/annihilation operators should 
be normalized so that 

[ak,al,] = 5kk', [a/c,afe/] = [al,al,] = 0 


The matrix Zij is integer and skew-symmetric, but oth¬ 
erwise arbitrary. Next, let 

M 

c' = '£v.,c,+x^ (13) 

i=i 

for some matrix V and some vector y. Then, [C[,C'j\ = 
2TTiZ[j where Z' = VZV^. The second step of the recipe 
is to find a matrix V with integer entries and determinant 
±1, such that Z' takes the simple form 






^ di 

0 . 


/07 

-V 



0 

(^2 ■ 

. 0 

V 

0 / 

0 

v = 




\0 

0 

0 

1 

to 










0 . 

• dij 


Here I is some integer with 0 < / < M/2 and 0/ denotes 
an J X / matrix of zeros. In mathematical language, 
V is an integer change of basis that puts Z into skew- 
normal form. It is known that such a change of basis 
always exists, though it is not unique.^® After finding an 
appropriate V, the offset y should then be chosen so that 

Xj = TT • ^ VijVikZjk (mod 2tt) (15) 

j<k 

The reason for choosing y in this way is that it ensures 
that = IV') for any |V') S HeS, as can be easily 

seen from the Campbell-Baker-Hausdorff formula. 

Once we perform these two steps, we can obtain the 
complete energy spectrum of iLgff with the help of a few 
results that we prove in appendix B.®® Our first result is 
that iLeff can always be written in the form 

K 

Lleff — E -f F ((727+1,..., (16) 

k=l 

where F is some (a priori unknown) quadratic function. 
Our second result (which is really just an observation) is 
that the following operators all commute with each other: 


We are now ready to discuss the second step of the 
recipe. This step involves searching for linear combi¬ 
nations of {Cl, ...,Cm} that have simple commutation 
relations with one another. The idea behind this step 
is that we ultimately need to construct a complete set 
of quantum numbers for labeling the eigenstates of Hgg. 
Some of these quantum numbers will necessarily involve 
the Ci operators since these operators play a prominent 
role in the definition of the effective Hilbert space, Hee- 
However, the (7i’s are unwieldy because they have com¬ 
plicated commutation relations with one another. Thus, 
it is natural to look for linear combinations of Ci’s that 
have simpler commutation relations. 

With this motivation in mind, let Zij be the M x M 
matrix defined by 


(12) 


,..., ,..., , C''7+i,..., C'm} 

(17) 

Furthermore, these operators commute with the occupa¬ 
tion number operators {a|ai,..., Therefore, we 

can simultaneously diagonalize (17) along with {a\ak}- 
We denote the simultaneous eigenstates by 


|6»i, ifi, X7+7, ...,x'jv7_7,ni, ...,nK) 

or, in more abbreviated form, \0,ip,x',n). Here the dif¬ 
ferent quantum numbers are defined by 

e<M|0,y,,a;',n) = e*®’|0, y), x', n), i = 1,..., J 
e'-^'\9,ip,x' ,n) = e'''^^-‘\6,ip,x' ,n), i = I 1,...,2/ 
C{\9,(p,x',n) = 2TTxl_j\9,(p,x',n), i = 2I+ 
alak\9,(p,x',n) =nk\9,(p,x',n), k = l,...,K (18) 
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where 0 < < 27r, while x[ is real valued and Uk 

ranges over non-negative integers. 

By construction the |0, a:', n) states form a complete 

basis for the Hilbert space T-L. Our third result is that 
a subset of these states form a complete basis for the 
effective Hilbert space Hee- This subset consists of all 
\0, ip, x', n) for which 

1.0 = (27rai/(ii,..., 2Trai/di) with ai = 0, 1 ,..., di — 1. 

2 . ¥>=( 0 , 0 ,..., 0 ). 

3. ...,x'm_j) = {qi, ...,qM-2i) for some integers 

qi- 

We will denote this subset of eigenstates by {|q:, q, n)}. 

Putting this together, we can see from equations (16) 
and (18) that the \cx,q,n) are eigenstates of HeS, with 
eigenvalues 

K 

E = '^nkEu+F{2TTqi,...,2'KqM-2i) (19) 

fc=i 

We therefore have the full eigenspectrum of HeS — up 
to the determination of the function F. With a bit more 
work, one can go further and compute the function F 
(see appendix B 3) but we will not discuss this issue here 
because in many cases of interest it is more convenient 
to find F using problem-specific approaches. 

To see examples of this diagonalization procedure, we 
refer the reader to section HI. As for the general deriva¬ 
tion of this procedure, see appendix B. 


Thus, the above formula predicts that the degeneracy 
for this system is D = y^det{Z) = d, which is consistent 
with our previous discussion. 


D. Finite U corrections 

We now discuss our last major result. To understand 
this result, note that while Hgff gives the exact low en¬ 
ergy spectrum of H in the infinite U limit, it only gives 
approximate results when U is large but finite. Thus, to 
complete our picture we need to understand what types 
of corrections we need to add to Heff to obtain an exact 
effective theory in the finite U case. 

It is instructive to start with a simple example: F[ = 

^ ^ - U cos(27ra;). As we discussed in section HA, 

the low energy effective Hamiltonian in the infinite U 
limit is iJeft = while the low energy Hilbert space 
HeS is spanned by position eigenstates {\q)} where q is 
an integer. 

Let us consider this example in the case where U is 
large but finite. In this case, we expect that there is some 
small amplitude for the system to tunnel from one cosine 
minima x = q to another minima, x = q — n. Clearly 
we need to add correction terms to 'HeS that describe 
these tunneling processes. But what are these correction 
terms? It is not hard to see that the most general possible 
correction terms can be parameterized as 


C. Degeneracy 

One implication of Eq. 19 which is worth mentioning 
is that the energy E is independent of the quantum num¬ 
bers ai, Since ai ranges from 0 < Oi < di — 1, it 

follows that every eigenvalue of Heff bas a degeneracy of 
(at least) 

I 

D = \{d, (20) 

i=l 

In the special case where Zij is non-degenerate (i.e. the 
case where M = 21), this degeneracy can be conveniently 
written as 

D = \/det(Z) (21) 


^ e“P-e„(x) (22) 

n——oo 

where en{x) is some unknown function which also de¬ 
pends on U. Physically, each term describes a tun¬ 
neling process \q) —)■ \q — n) since e^^P\q) = \q — n). The 
coefficient e„(a;) describes the amplitude for this process, 
which may depend on q in general. (The one exception is 
the n = 0 term, which does not describe tunneling at all, 
but rather describes corrections to the onsite energies for 
each minima). 

Having developed our intuition with this example, we 
are now ready to describe our general result. Specifically, 
in the general case we show that the finite U corrections 
can be written in the form 


since 

I 

det(Z) = det(Z') = 

i=l 


For an example of this degeneracy formula, consider 
the Hamiltonian (9) discussed in section HA. In this case. 
Cl = dp while C 2 = 2ttx so 





(23) 


with the sum running over M component integer vectors 
m = (mi, ...,mM)- Here, the are unknown functions 
of {oi,..., Ofc, a|,..., a\, C 21 + 1 ', Cm} which also depend 
on U. We give some examples of these results in section 
HI. For a derivation of the finite U corrections, see ap¬ 
pendix C. 
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E. Splitting of ground state degeneracy 


One application of Eq. 23 is to determining how the 
ground state degeneracy of i?eff splits at finite U. Indeed, 
according to standard perturbation theory, we can find 
the splitting of the ground state degeneracy by projecting 
the finite U corrections onto the ground state subspace 
and then diagonalizing the resulting D x D matrix. The 
details of this diagonalization problem are system depen¬ 
dent, so we cannot say much about it in general. How¬ 
ever, we would like to mention a related result that is 
useful in this context. This result applies to any system 
in which the commutator matrix Zij is non-degenerate. 
Before stating the result, we first need to define some 
notation: let Ti, ...jTm be operators defined by 

= (24) 

j 


Note that, by construction, the T^ operators obey the 
commutation relations 


[Q,r,]=27rz%, [afc,r,] = [4,r,]=0 (25) 


With this notation, our result is that 


a e 


a = Mr 


(a e 


E m 

j=i’ 


|a) (26) 


where jo;), la') are ground states and Um is some un¬ 
known proportionality constant. This result is useful be¬ 
cause it is relatively easy to compute the matrix elements 
of e ^ ^; hence the above relation allows us to com¬ 

pute the matrix elements of the finite U corrections (up 
to the constants Um) without much work. We derive this 
result in appendix C. 


III. EXAMPLES 

In this section, we illustrate our formalism with some 
concrete examples. These examples involve a class of 
two dimensional electron systems in which the spin-up 
and spin-down electrons form v = 1/k Laughlin states 
with opposite chiralities^^. These states are known as 
“fractional quantum spin Hall insulators.” We will be 
primarily interested in the edges of fractional quantum 
spin Hall (FQSH) insulators^^. Since the edge of the 
Laughlin state can be modeled as a single chiral Luttinger 
liquid, the edge of a FQSH insulator consists of two chiral 
Luttinger liquids with opposite chiralities — one for each 
spin direction (Fig. I). 

The examples that follow will explore the physics of the 
FQSH edge in the presence of impurity-induced scatter¬ 
ing. More specifically, in the first example, we consider 
a FQSH edge with a single magnetic impurity; in the 
second example we consider a FQSH edge with multiple 
magnetic impurities; in the last example we consider a 
FQSH edge with alternating magnetic and superconduct¬ 
ing impurities. In all cases, we study the impurities in the 



FIG. 1. The fractional quantum spin Hall edge consists of 
two counter-propagating chiral Luttinger Liquids — one for 
each spin direction (t,4) 


infinite scattering limit, which corresponds to f7 —>■ oo in 
(1). Then, in the last subsection we discuss how our re¬ 
sults change when the scattering strength U is large but 
finite. 

We emphasize that the main purpose of these examples 
is to illustrate our formalism rather than to derive novel 
results. In particular, many of our findings regarding 
these examples are known previously in the literature in 
some form. 

Of all the examples, the last one, involving magnetic 
and superconducting impurities, is perhaps most interest¬ 
ing: we find that this system has a ground state degen¬ 
eracy that grows exponentially with the number of im¬ 
purities. This ground state degeneracy is closely related 
to the previously known topological degeneracy that ap¬ 
pears when a FQSH edge is proximity coupled to alter¬ 
nating ferromagnetic and superconducting strips^^^^^. 

Before proceeding, we need to explain what we mean 
by “magnetic impurities” and “superconducting impuri¬ 
ties.” At a formal level, a magnetic impurity is a local¬ 
ized object that scatters spin-up electrons into spin-down 
electrons. Likewise a superconducting impurity is a lo¬ 
calized object that scatters spin-up electrons into spin- 
down holes. More physically, a magnetic impurity can be 
realized by placing the tip of a ferromagnetic needle in 
proximity to the edge while a superconducting impurity 
can be realized by placing the tip of a superconducting 
needle in proximity to the edge. 


A. Review of edge theory for clean system 

As discussed above, the edge theory for the v = Ijk 
fractional quantum spin Hall state consists of two chi¬ 
ral Luttinger liquids with opposite chiralities — one for 
each spin direction (Fig. 1). The purpose of this section 
is to review the Hamiltonian formulation of this edge 
theory.^’^’^^ More specifically we will discuss the edge 
theory for a disk geometry where the circumference of 
the disk has length L. Since we will work in a Hamil¬ 
tonian formulation, in order to define the edge theory 
we need to specify the Hamiltonian, the set of physical 
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observables, and the canonical commutation relations. 

We begin with the set of physical observables. 
The basic physical observables in the edge theory 
are a collection of operators {dy(l)^{y),dy(j)i{y)} along 
with two additional operators (p^{yo), ipliyo) where 
yo is an arbitrary, but fixed, point on the bound¬ 
ary of the disk. The {dy(j>^{y),dy(j)i{y),(j)^{yo),(l)i{yo)} 
operators can be thought of as the fundamental 
phase space operators in this system, i.e. the ana¬ 
logues of the {xi, ...,xn,Pi, ■■■,Pn} operators in sec¬ 
tion II. Like {xi, ...,xn,Pi, ...,pn}, all other physical 
observables can be written as functions/functionals of 
{dy(j)^{y),dy(l>i{y),cl)^{yo),(l>l{yo)}. Two important ex¬ 
amples are the operators (j)^{y) and 4>i{y) which are de- 
fined by 

ry 

(j)a{y) = (jicriyo) + / dx(j)adx, cr=t,4. (27) 

dyo 


usual phase space operators {xi, ...,xn,Pi, ■■■,Pn}, the 
commutators of {dy^^{y),dy(j)i{y),(p^{yo),(j)i{yQ)} are c- 
numbers. More specifically, the basic commutation rela¬ 
tions are 

Stt? 

[da:(l}^{x), dy(j)^{y)] = —dj{x - y) 

O'TTO 

[d^(l3^{x), dy(j)i{y)\ = —-^dJix - y) 

O-TT? 

[</>T(2/o),9y<(>t(2/)] = -Y^(y~yo) 

Ott? 

Ihiyo), dy(l)i{y)] = — —S{y - yo) (28) 

with the other commutators vanishing: 

[4’t(.yo),dy(j)i{y)] = [<));(yo),9y(^t(y)] = 0 

[Myo),hiyo)] = [dxM^)^9y(i)i{y)] = o 


where the integral runs from yo to y in the clockwise 
direction. 

The physical meaning of these operators is as follows: 
the density of spin-up electrons at position y is given by 
Pt(?/) = while the density of spin-down electron is 

Piiv) = The total charge Q and total spin on 

the edge are given by Q = and = I/2((5t—Q^) 

with 

^ /•L /2 

Qa — J dy(f)ady, (T ="(, 4 - 

Finally, the spin-up and spin-down electron creation op¬ 
erators take the form 

In the above discussion, we ignored an important sub¬ 
tlety: 4>^{yo) and 4>i[yo) are actually compact degrees of 
freedom which are only defined modulo 27r/fc. In other 
words, strictly speaking, 4)^{yo) and (('^(yo) are not well- 
defined operators: only and are well- 

defined. (Of course the same also goes for (j)^{y) and 
4>i{y), in view of the above definition). Closely related 
to this fact, the conjugate “momenta” and are 
actually discrete degrees of freedom which can take only 
integer values. 

The compactness of 4’^{yo)i4‘i{yQ) and discreteness of 
Qt ’ Qi inconvenient for us since the machinery dis¬ 
cussed in section II is designed for systems in which all 
the phase space operators are real-valued, rather than 
systems in which some operators are angular valued and 
some are integer valued. To get around this issue, we 
will initially treat 4>^{yo) and ())^(yo) and the conjugate 
momenta Qtj Qi valued operators. We will then 

use a trick (described in the next section) to dynamically 
generate the compactness of 4'tiyo)^4‘liyo) as well as the 
discreteness of Q^,Qi. 

Let us now discuss the commutation relations for the 
{dy(j)^{y),dy(j>i{y),(j)^{yQ),(j>l{yo)} operators. Like the 


Using these basic commutation relations, together with 
the definition of (t>a{y) (27), one can derive the more gen¬ 
eral relations 

O-TT? 

[(p^ix), dycp^iy)] = —S{x-y) 

Ott? 

[(plix), dycpiiy)] = — —6{x - y) 
[(j)^{x),dy(j)i{y)] = 0 (29) 

as well as 

TTZ 

[hi^)^hiy)] = -^sgn{x,y) 

TTZ 

[(Plix),(l)i{y)] = -—sgii{x,y) 
[(j)^{x),(j)i{y)]=0 (30) 

where the sgn function is defined by sgn(a;,y) = -Tl if 
yo < X < y and sgn(a:,y) = —\ \i yo < y < x, with the 
ordering defined in the clockwise direction. The latter 
commutation relations (29) and (30) will be particularly 
useful to us in the sections that follow. 

Having defined the physical observables and their com¬ 
mutation relations, the last step is to define the Hamilto¬ 
nian for the edge theory. The Hamiltonian for a perfectly 
clean, homogeneous edge is 

Ho = ^ [{dxM^)f + {dxM^)f]dx (31) 

47r 

where v is the velocity of the edge modes. 

At this point, the edge theory is complete except for 
one missing element: we have not given an explicit defini¬ 
tion of the Hilbert space of the edge theory. There are two 
different (but equivalent) definitions that one can use. 
The first, more abstract, definition is that the Hilbert 
space is the unique irreducible representation of the oper¬ 
ators {dy(l)^,dy(l)i,(j)^{yo),4>i{yQ)} and the commutation 
relations (28). (This is akin to defining the Hilbert space 


cosine terms to our Hamiltonian: 




FIG. 2. (a) A magnetic impurity on a fractional quantum 
spin Hall edge causes spin-up electrons to backscatter into 
spin-down electrons, (b) In the infinite backscattering limit, 
the impurity effectively reconnects the edge modes. 


of the ID harmonic oscillator as the irreducible represen¬ 
tation of the Heisenberg algebra [x,p\ = i). The second 
definition, which is more concrete but also more compli¬ 
cated, is that the Hilbert space is spanned by the com¬ 
plete orthonormal basis {l?^, {''T'pf})} where the 

quantum numbers 9 ^, 9 ^ range over all integers^® while 
rip^, Upi range over all nonnegative integers for each value 
of p = 27r/L,47r/L,.... These basis states have a simple 
physical meaning: \q^,qi,{np-f},{npi}) corresponds to a 
state with charge and q^ on the two edge modes, and 
with Up-f and Hpi phonons with momentum p on the two 
edge modes. 


B. Example 1: Single magnetic impurity 

With this preparation, we now proceed to study a 
fractional quantum spin Hall edge with a single mag¬ 
netic impurity in a disk geometry of circumference L 
(Fig. 2a). We assume that the impurity, which is lo¬ 
cated at x = 0 , generates a backscattering term of the 
form ^('0|(O)^^(O) -I- h.c). Thus, in the bosonized lan¬ 
guage, the system with an impurity is described by the 
Hamiltonian 

H = Ho-UcosiC), C = fc((/.t(0) + </>;(0)) (32) 

where Hq is defined in Eq. 31. Here, we temporarily 
ignore the question of how we regularize the cosine term; 
we will come back to this point below. 

Our goal is to find the low energy spectrum of H in the 
strong backscattering limit, U —>■ 00 . We will accomplish 
this using the results from section H. Note that, in using 
these results, we implicitly assume that our formalism 
applies to systems with infinite dimensional phase spaces, 
even though we only derived it in the finite dimensional 
case. 

First we describe a trick for correctly accounting for 
the compactness of (/>t( 2 /o)) (j/o) and the quantization 

of The idea is simple: we initially treat these 

variables as if they are real valued, and then we intro¬ 
duce compactness dynamically by adding two additional 


H = Hq — U cos(C') — U cos(27r(5t) ~ U cos(27r(3^) (33) 

These additional cosine terms effectively force Q-f and 
Ql to be quantized at low energies, thereby generating 
the compactness that we seek.^® We will include all three 
cosine terms in our subsequent analysis. 

The next step is to calculate the low energy effective 
Hamiltonian Ffgff and low energy Hilbert space 'HeS- In¬ 
stead of working out the expressions in Eqs. 2, 5, we 
will skip this computation and proceed directly to find¬ 
ing creation and annihilation operators for HeS using Eq. 
11. (This approach works because equation (11) does not 
require us to find the explicit form of Hgs). 

According to Eq. 11, we can find the creation and anni¬ 
hilation operators for Hgs by finding all operators a such 
that (1) a is a linear combination of our fundamental 
phase space operators {dy4i^,dy(j)i,4‘^iyo),4>iiyo)} and 
(2) a obeys 

[a, Hq] = Ea + Hq] + A^[(5^, Hq] + A^[(3^, Hq] 

[a, C] = [a, Q^] = [a, Q;] = 0 (34) 

for some scalars A, A, A^, A^ with E y^O. 

To proceed further, we note that the constraint 
[a,Q^] = [a, Q;] = 0, implies that <(>t(yo), <?i’;(2/o) cannot 
appear in the expression for a. Hence, a can be written 
in the general form 

rL/2 

a= [f^{y)dy(l)^{y) + fi{y)dycl)^{y)]dy (35) 
J-L/2 

Substituting this expression into the first line of Eq. 34 
we obtain the differential equations 

-ivf^iy) = Ef^{y) + Xkiv5{y) 
ivf[{y) = Ef^iy) - XkivSiy) 

(The A^,A^ terms drop out of these equations since 
commute with Hq). These differential equations 
can be solved straightforwardly. The most general solu¬ 
tion takes the form 

/t(2/)=e*™[Ai0 (-j/)+A20 (y)] 

h{y) = e-^Py[B,Q{-y) + i?20(2/)] (36) 

where p = E/v^ and 

A2 = Ai- Afc, B2 = Bi- Xk (37) 


Here 0 is the Heaviside step function defined as 


0(x) 


0 -L/2 <x<0 
1 0 < X < L/2 


(Note that the above expressions (36) for /^, do not 
obey periodic boundary conditions at x = ±L/2; we will 
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not impose these boundary conditions until later in our 
calculation). Eliminating A from (37) we see that 

A 2 — Ai = B 2 — Bi (38) 

We still have to impose one more condition on a, 
namely [a, C] = 0. This condition leads to a second 
constraint on Ai, A 2 , Bi^ B 2 , but the derivation of this 
constraint is somewhat subtle. The problem is that if we 
simply substitute (35) into [a, C] = 0, we find 

/t(0) = 4(0) (39) 

It is unclear how to translate this relation into one in¬ 
volving Ai, A 2 , Bi, B 2 since 4,/t discontinuous at 
X = 0 and hence 4(0), 4(0) are ill-defined. The origin 
of this puzzle is that the cosine term in Eq. 32 contains 
short-distance singularities and hence is not well-defined. 
To resolve this issue we regularize the argument of the 
cosine term, replacing C = fc(4(0) + 4(0)) with 

/•V2 

C —>■ / k((:j)^(x) + (j)i{x))d{x)dx (40) 

J-L/2 


where p is quantized as p = ±tt/L, ±2Tr/L ,... and Ep = 
vp. (Note that p = 0 does not correspond to a legitimate 
creation/annihilation operator according to the definition 
given above, since we require E ^ 0). 

Following the conventions from section IIB, we will re¬ 
fer to the operators with Ep > 0 — or equivalently p > 0 
— as “annihilation operators” and the other operators as 
“creation operators.” Also, we will choose the normal¬ 
ization constant Ai so that [ap,ap,] — Spp/ for p,p' > 0. 
This gives the expression 

I k 

y_^^^(e-44 + e--44)e(-4 

(44) 

The next step is to compute the commutator matrix 
Zij. In the case at hand, we have three cosine terms 
{cos(C'i),cos(C'2),cos(C'3)}, where 

Cl =C, C2= 2nQ^, C 3 = 2ttQ^ 


where <5 is a narrowly peaked function with f 5{x)dx = I. 
Here, we can think of 5 as an approximation to a delta 
function. Note that <5 effectively introduces a short- 
distance cutoff and thus makes the cosine term non¬ 
singular. After making this replacement, it is straight¬ 
forward to repeat the above analysis and solve the differ¬ 
ential equations for 4; 4- appendix G, we work out 
this exercise, and we find that with this regularization, 
the condition [a, C] = 0 leads to the constraint 

Aidi = ?i±Ei ,41) 

Combining our two constraints on Ai, Hi, A2, B 2 (38,41), 
we obtain the relations 

Ai=Bi, A2=B2 (42) 

So far we have not imposed any restriction on the mo¬ 
mentum p. The momentum constraints come from the 
periodic boundary conditions on 4i 4- 

M-L/2) = 4(L/2), hi-L/2) = 4(L/2) 

Using the explicit form of 4)4) these boundary condi¬ 
tions give 

^ BiEP^C = B2e-^P^C 

from which we deduce 


Therefore Z^j is given by 



1 

0 

0 



To proceed further we need to find an appropriate 

change of variables of the form C' = '^ijCj + Xi- 

Here, V should an integer matrix with determinant ±1 
with the property that Zb = 2¥l[^()^i] skew- 

normal form, while x should be a real vector satisfying 
Eq. 15. It is easy to see that the following change of 
variables does the job: 

C[= Cl, C'2 = -27rQ^, C'^ = 27rQ^ + 2^Q; 

Indeed, for this change of variables, 

/O -1 0 
Z' = 1 0 0 

\0 0 0 

We can see that this is in the canonical skew normal form 
shown in Eq. 14, with the parameters M = 3, / = I, 
di = 1 . 

We are now in a position to write down the low energy 
effective Hamiltonian i7eff: according to Eq. 16, HeS 
must take the form 


g2ipL ^ ^ Aie-*P^ (43) 

Putting this all together, we see that the most general 
possible creation/annihilation operator for EleS is given 
by 

rL/2 

ap = Ai / (e*P«44 + e-*P^44)0(-p) 

J-L/2 

+e-^pL{EPydy^^ + e-*™44)0(y)dy 


HeS = ^ vpalap + F ■ (C'^f (45) 

p>0 

where F is some (as yet unknown) constant. To deter¬ 
mine the constant F, we make two observations. First, 
we note that the first term in Eq. 45 can be rewritten 
as J2p^o ^(“-pOp)- Second, we note that Cg = 27rQ is 
proportional to Op^g- Given these observations, it is nat¬ 
ural to interpret the F(C^)^ term as the missing p = 0 
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term in the sum. This suggests that we can fix the co¬ 
efficient F using continuity in the p —> 0 limit. To this 
end, we observe that 

We conclude that F = Substituting this into (45), 
we derive 

Hes = '^p4>^p + ^ • ( C ' a )^ ( 46 ) 

p>0 




where the sum runs over p = njL, 27r/L,.... 

In addition to the effective Hamiltonian, we also need 
to discuss the effective Hilbert space HeS in which this 
Hamiltonian is defined. According to the results of sec¬ 
tion HB, the effective Hilbert space T-LeS is spanned by 
states {I?, {up})} where \q,{np}) is the unique simulta¬ 
neous eigenstate of the form 

F’^'^\q,{np}) = |g,{np}), 

6**^2 |g, {rip}) = |g,{np}), 

Cslq, {np}) = 2TTq\q, {rip}), 
alap\q, {rip}) = np\q, {rip}) 

Here Up runs over non-negative integers, while q runs 
over all integers. Note that we do not need to label the 
{b, {rip})} basis states with a quantum numbers since 
di = 1 so there is no degeneracy. 

Having derived the effective theory, all that remains 
is to diagonalize it. Fortunately we can accomplish this 
without any extra work: from (46) it is clear that the 
{|g, {rip})} basis states are also eigenstates of iJeS with 
energies given by 

E = Y^vpnpF — -q (47) 

p>0 

We are now finished: the above equation gives the com¬ 
plete energy spectrum of ideff, and thus the complete low 
energy spectrum of FI in the limit [/ —> oo. 

To understand the physical interpretation of this en¬ 
ergy spectrum, we can think of rip as describing the num¬ 
ber of phonon excitations with momentum p, while q de¬ 
scribes the total charge on the edge. With these identifi¬ 
cations, the first term in (47) describes the total energy 
of the phonon excitations — which are linearly dispers¬ 
ing with velocity v — while the second term describes 
the charging/capacitative energy of the edge. 

It is interesting that at low energies, our system has 
only one branch of phonon modes and one charge de¬ 
gree of freedom, while the clean edge theory (31) has two 
branches of phonon modes and two charge degrees of free¬ 
dom — one for each spin direction. The explanation for 
this discrepancy can be seen in Fig. 2b: in the infinite U 
limit, the impurity induces perfect backscattering which 
effectively reconnects the edges to form a single chiral 
edge of length 2L. 


FIG. 3. (a) A collection of N magnetic impurities on a frac¬ 
tional quantum spin Hall edge. The impurities are located at 
positions x\,...,xn- (b) In the infinite backscattering limit, 
the impurities effectively reconnect the edge modes, breaking 
the edge into N disconnected components. 


C. Example 2: Multiple magnetic impurities 


We now consider a fractional quantum spin Hall edge 
in a disk geometry with N magnetic impurities located at 
positions cci,..., a:AT (Fig. 3a). Modeling the impurities in 
the same way as in the previous section, the Hamiltonian 
is 

N 

H = cos(C'j), Ci = 

(48) 

where F[q is defined in Eq. 31. 

As in the single impurity case, our goal is to understand 
the low energy physics of H in the limit U ^ oo. We 
can accomplish this using the same approach as before. 
The first step is to take account of the compactness of 
the discrete nature of Q^,Qi by adding two 
additional cosine terms to our Hamiltonian: 

N 

F[ = Hq — U ^ cos(Ci) — U cos{2ttQ^) — U cos{2ttQi) 

i=l 

Next, we find the creation and annihilation operators 
for Hes using Eq. 11. That is, we search for all operators 
a such that (1) a is a linear combination of our fundamen¬ 
tal phase space operators {dy4>^,dy4>i,(j)^{yo),4‘i{yo)} 
and (2) a obeys 


N 

[a, Hq] = Ea + '^2, Hq] , Hq] + , F[q\ 

i=l 


[a, Cj] = [a, Q^] = [a, Q;] = 0 


(49) 


for some E, A^, A^, A^ with E ^ 0. 

Given that [a, Q^] = [a, Qi] = 0, we know that 
4>t{yo)j4>liyo) cannot appear in the expression for a. 
Hence, a can be written in the general form 


a = 



[hiy)9yH(.y) + hiy)dyMy)]dy 
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Substituting these expressions into the first line of Eq. 49, 
we obtain 


N 

= Ef^{y) + kiv ^ \j5{y - xj) 
i=i 

N 

ivfliy) = Efi{y) - kiv'^Xj6{y - Xj) 

1=1 

Solving these differential equations gives 

N 

ftiy) = <y< Xj)] 

1=1 

N 

hiy) = Bje~"pye{xj -1 <y< Xj)] 

1=1 

where p = E/v and where 

Aj+i = Aj - \jke-^P^^ , Bj+i = Bj - XjkEP^^ 

Here, Q{a < y < b) is defined to take the value 1 if y 
is in the interval [o, b] and 0 otherwise. Also, we use a 
notation in which xq is identified with xn. Eliminating 
Xj, we derive 

(A,+i - (50) 

We still have to impose the condition [a, Cj] = 0, which 
gives an additional set of constraints on {Aj, Bj}. As in 
the single impurity case, we regularize the cosine terms 
to derive these constraints. That is, we replace C = 
k{ip^{xj) + 'Piixj)) with 

/•V 2 

C —i / k{(j)^{x) + 4)i{x))S{x — Xj)dx (51) 

J-L/2 

where (5 is a narrowly peaked function with J 5{x)dx = 1 , 
i.e., an approximation to a delta function. With this 
regularization, it is not hard to show that [a, Cj] = 0 
gives the constraint 

i(A, + = \{B, + (52) 

Combining (50), (52) we derive 

AjCP^^ = , Aj+iEP^^ = B^+ie-^P^^ (53) 

Our task is now to find all {Aj, Bj,p} that satisfy 
(53). For simplicity, we will specialize to the case where 
the impurities are uniformly spaced with spacing s, i.e. 
Xj+i —Xj = s = L/N for all j. In this case, equation (53) 
implies that = 1 , so that p is quantized in integer 
multiples of tt/s. For any such p, (53) has N linearly 
independent solutions of the form 


with m = 1, Putting this all together, we see that 

the most general possible creation/annihilation operator 
for Hgg is given by 

I k 

• Q{Xm-l <y < Xm)\dy 


with Epm = vp. Here the index m runs over to = 1,..., A 
while p takes values ±7r/s, ±27r/s,.... (As in the single 
impurity case, p = 0 does not correspond to a legitimate 
creation/annihilation operator, since we require that E ^ 

0 ). ^ 

Following the conventions from section HB, we will 
refer to the operators with E > {) — or equivalently p > 0 
— as “annihilation operators” and the other operators as 
“creation operators.” Note that we have normalized the 
a operators so that [apmjOp,^,] = Spp'Smm' for p,p' > 0. 

The next step is to compute the commutator matrix 
Zij = ^ [Ci,Cj\. Let us denote the A+2 cosine terms as 
{cos(C'i), ...,cos(Cjv+ 2 )} where Cn+i = 27rQ^, Cn +2 = 
iTrQj. Using (30) we find that Zij takes the form 


/ 0 ••• 0 1 -i\ 


Zij — 


0 ••• 0 1-1 
-1 ••• -10 0 
V 1 ••• 1 0 0 / 


To proceed further we need to find an appropriate 
change of variables of the form C' = + Xi- 

Here, V should be chosen so that Zk = 2 ^[C'',C''] is in 
skew-normal form, while x should be chosen so that it 
obeys Eq. 15. It is easy to see that the following change 
of variables does the job: 


= Cl, C 2 = -27rQ^, C 3 = 27rQ^ 2ttQi 
C!^ = Cm-2 ~ Cm-3) TO = 4, ..., N 2 


0 

Vo 0 0 ••• 0 / 

We can see that this is in the canonical skew-normal form 
shown in Eq. 14, with the parameters M = N + 2, I = 1, 
di = 1 . 

We are now in a position to write down the low energy 
effective Hamiltonian HeS'. according to Eq. 16, Hgg 
must take the form 

N 

HeS='^Yl '^P^lmapm + ..., C^v+s) (54) 

m—1 p>0 


Indeed, it is easy to check that 


/O -1 0 • 
' 1 0 0 • 
0 0 0 - 


Bj = Aj = 0 for j m 
Bj = Aje^^P^^ ^ 0 for j = to 
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where the sum runs over p = tt/s, 27r/s,... and where F 
is some quadratic function of N variables. To determine 
F, we first need to work out more concrete expressions 
for C'j^. The m = 3 case is simple: C'^ = 2TrQ. On the 
other hand, for m = 4,..., + 2, we have 

Cm = 

m—3) 

l‘Xm-2 

= k {dyc/)^ + dy(l)i)dy 

J Xm-3 


where the second line follows from the definition of 0^, (f>i 
(27) along with the assumption that the impurities are 
arranged in the order yo < xi < ... < xn in the clockwise 
direction. 

With these expressions we can now find F. We use 
the same trick as in the single impurity case: we note 
that the first term in Eq. 54 can be rewritten as as 
Sm J2p^o ^(a-pmOpm), and we observe that 


r vIpI 
hm -—a-pma 
p-s-O Z 


pm 


V 

Snks 


icL+2r 


for m = 2,..., N, while 


v\p\ 


V 

Sirks 


Ar+2 

• (fcC' - ^ Cm'f 

m—4 


Assuming that F reproduces the missing p = 0 piece of 
the first term in Eq. 54, we deduce that 

Af+2 

m—4: 

iV+2 

+ ^ -1: (») 

m—4 

In addition to the effective Hamiltonian, we also need 
to discuss the effective Hilbert space HeS- Applying the 
results of section HB, we see that T-LeS is spanned by 
states {|<7, {upm})} where \q, {upm}) is the unique simul¬ 
taneous eigenstate of the form 

e*C|q^{„p^}) = \q,{npm}): 
e ^1?, {Upm}) = |9;{^pm})) 

Cj \ q: {n^pm }) — ‘^'^qi—2 I q7 {^pm}) 5 ^ • 7 ^ F 2 

apm^pml^i {^pm}) — ^pml^i {^pm}) 

Here Upm runs over non-negative integers, while q is an 
N component vector, q — (gi, ...,qN) where each compo¬ 
nent qi runs over all integers. As in the single impurity 
case, we do not need to label the {\q, {upm})} basis states 
with a quantum numbers since di = 1 and thus there is 
no degeneracy. 

Now that we have derived the effective theory, all that 
remains is to diagonalize it. To do this, we note that the 


{|q, {^pm})} basis states are also eigenstates of iJeS with 
energies given by 

AT AT 

c = ^ ^vpup^ + ^ H (56) 

m—1 p>0 m—2 

+ - q2 - - - qN) 

The above equation gives the complete low energy spec¬ 
trum of F[ in the limit C/ —> oo. 

Let us now discuss the physical interpretation of these 
results. As in the single impurity case, when f7 —>■ oo 
the impurities generate perfect backscattering, effectively 
reconnecting the edge modes. The result, as shown in 
Eig. 3b, is the formation of N disconnected chiral modes 
living in the N intervals, [a::Ar,a:i], [xi, 0:2],..., [xAr_i,XAr]. 

With this picture in mind, the Upm quantum numbers 
have a natural interpretation as the number of phonon 
excitations with momentum p on the mth disconnected 
component of the edge. Likewise, if we examine the def¬ 
inition of qm: we can see that qm/k is equal to the total 
charge in the mth component of the edge, i.e. the total 
charge in the interval [xm-i,Xm], for m = 2,...,N. On 
the other hand, the quantum number qi is slightly differ¬ 
ent: qi is equal to the total charge on the entire boundary 
of the disk [—L/2, L/2]. Note that since qm is quantized 
to be an integer for all m = 1, ...,N, it follows that the 
charge in each interval [xm-i,Xm] is quantized in integer 
multiples of 1/k while the total charge on the whole edge 
is quantized as an integer. These quantization laws are 
physically sensible: indeed, the fractional quantum spin 
Hall state supports quasiparticle excitations with charge 
1/fc, so it makes sense that disconnected components of 
the edge can carry such charge, but at the same time we 
also know that the total charge on the boundary must be 
an integer. 

Putting this all together, we see that the first term in 
(56) can be interpreted as the energy of the phonon exci¬ 
tations, summed over all momenta and all disconnected 
components of the edge. Similarly the second term can 
be interpreted as the charging energy of the disconnected 
components labeled by m = 2,..., N, while the third term 
can be interpreted as the charging energy of the first com¬ 
ponent labeled by m = 1. 

So far in this section we have considered magnetic im¬ 
purities which backscatter spin-up electrons into spin- 
down electrons. These impurities explicitly break time 
reversal symmetry. However, one can also consider non¬ 
magnetic impurities which preserve time reversal symme¬ 
try and backscatter pairs of spin-up electrons into pairs 
of spin-down electrons. When the scattering strength U 
is sufficiently strong these impurities can cause a spon¬ 
taneous breaking of time reversal symmetry, leading to 
a two-fold degenerate ground state.This physics 
can also be captured by an appropriate toy model and 
we provide an example in Appendix H. 
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(a) (b) 

FIG. 4. (a) A collection of 2N alternating magnetic 

and superconducting impurities on a fractional quantum spin 
Hall edge. The magnetic impurities are located at posi¬ 
tions Xi,X 3 , ....,X 2 N-i while the superconducting impurities 
are located at positions X 2 ,X 4 ,, ...,X 2 n- The magnetic impu¬ 
rities scatter spin-up electrons into spin-down electrons while 
the superconducting impurities scatter spin-up electrons into 
spin-down holes, (b) In the infinite U limit, the impurities 
effectively reconnect the edge modes, breaking the edge into 
2N disconnected components. 



D. Example 3: Multiple magnetic and 
superconducting impurities 


and (2) a obeys 


2N 


[o, Hq] — Ea + WCi, Hq] , Hq] + \Q ^, 


i=l 


[a, Cj] = [a, Q^] = [a, Q^] = 0 


(60) 


for some E, with E ^ 0. 

As before, since [a, = [a, Qi\ = 0, it follows that 

4>^(yo),4‘i{yQ) cannot appear in the expression for a. 
Hence, a can be written in the general form 


a = 



[h{y)dyh{y) + fi{y)dvh{y)¥y 


Substituting this expression into the first line of Eq. 60, 
we obtain 

-ivf\{y) = Ef^ + kiv'^Xj6{y - xj) 

3 

ivf[{y) = Ef^ -kivY^{-iy+^XjS{y- Xj) 

3 

Solving the above first order differential equation we get 


We now consider a fractional quantum spin Hall edge 
in a disk geometry of circumference L with ‘IN alter¬ 
nating magnetic and superconducting impurities. We 
take the magnetic impurities to be located at positions 
xi, xa,..., a;2Ar-i while the superconducting impurities 
are located at positions x^.x/^, ...,X 2 n (Fig. 4a). We as¬ 
sume that the magnetic impurities generate a backscat- 
tering term of the form ^(^/^|(0)^/>^(0) -I- /i.c), while the 
superconducting impurities generate a pairing term of 
the form ^+ h.c). The Hamiltonian is then 

2N 

H = Ho- U'^cos{Cy (57) 

Ci = k{(l)^{xi) + (-1)*+V4(a:*)) (58) 

where Hq is defined in Eq. 31. 

As in the previous cases, our goal is to understand 
the low energy physics of H in the limit f7 —>■ oo. As 
before, we take account of the compactness of 4>^, 4>^ and 
the discrete nature of Q^,Q^ by adding two additional 
cosine terms to our Hamiltonian: 

2N 

H = Hq — U ^ cos(C'i) — U cos{2ttQi^) — U cos{2ttQi) 

i=l 

(59) 

Next, we find the creation and annihilation operators 
for Heff using Eq. 11. That is, we search for all operators 
a such that (1) a is a linear combination of our fundamen- 
tal phase space operators (/>^(2/o), 0;(yo)} 


N 

ftiy) = 

1=1 

N 

hiy) = ^Bje-'^PyQ{xj-i <y< Xj)] 

1=1 

where p = E/v and where 

Aj+i = Aj - Xjke-^P^EBj+i = Bj - {-ly+^XjkEP^^ 

Here, Q{a < y < b) is defined to take the value 1 if y 
is in the interval [a, b] and 0 otherwise. Also, we use a 
notation in which xq is identified with xn- Eliminating 
Xj, we derive 

(A,+1 - Ayyp^^ = (-1)^+1 (B.+i - Bye-^p^^ (61) 

We still have to impose the requirement [a,Cj\ = 0 
and derive the corresponding constraint on {Aj, Bj}. As 
in the previous cases, the correct way to do this is to 
regularize the cosine terms, replacing 

Cj—>■ / k{(j)j{x) + {—iy~^^(j>i{x))S{x — Xj)dx (62) 

J-L/2 

where (5 is a narrowly peaked function with J 5{x)dx = 1, 
i.e., an approximation to a delta function. With this 
regularization, it is not hard to show that [a, Cj] = 0 
gives the constraint 

^(A, + = \{-iy+yB, + B,+ye-^P^^ (63) 
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Combining (61), (63) we derive 

= {-ly+^Bje-^P^A 

= {-ly+^Bj+ie-^P^^ (64) 

Our task is now to find all {Aj, Bj,p} that satisfy 
(64). For simplicity, we will specialize to the case where 
the impurities are uniformly spaced with spacing s, i.e. 
Xj+i — Xj = s = L/N for all j. In this case, equation 
(53) implies that = —1, so that p is quantized in 
half-odd-integer multiples of tt/s. For any such p, (53) 
has N linearly independent solutions of the form 


Bj = Aj = {) for j y m 

Bj = {-ly+^AjC^^P^^ yO for j = m 


with m = 1,..., 2N. 

Putting this all together, we see that the most general 
possible creation/annihilation operator for H^g is given 

by 


^pm 







■ Q{Xjn-l <y < Xm)]dy ■ 



with Epm = vp. Here the index m runs over 
m = 1,...,2N while p takes values ± 7 r/ 2 s, ± 37 r/ 2 s,.... 
Note that we have normalized the a operators so that 

[^prm P-jP ^ 0- 

The next step is to compute the commutator ma¬ 
trix Zij = ^[Ci,Cj]. Let us denote the 2N + 2 co¬ 
sine terms as {cos(Ci),..., cos(C' 2 Ar+ 2 )} where C 2 N +1 = 
2TTQj, C 2 N +2 = 27rQ^. Using the commutation relations 
(30), we find 


/ 

0 

k 

0 

k ■■ 

• 0 

k 

1 



-k 

0 

k 

0 • • 

• k 

0 

1 

1 


0 

-k 

0 

k ■■ 

■ 0 

k 

1 

-1 


-k 

0 

-k 

0 • • 

• k 

0 

1 

1 


0 

-k 

0 

-k ■■ 

■ 0 

k 

1 

-1 


-k 

0 

-k 

0 • • 

• -k 

0 

1 

1 


-1 

-1 

-1 

-1 •• 

■ -1 

-1 

0 

0 

V 

1 

-1 

1 

-1 •• 

• 1 

-1 

0 

0/ 


To proceed further we need to find an appropriate 
change of variables of the form C' = + Xi- 

Here, V should be chosen so that Zy = ^[C[,Cj] is in 
skew-normal form, while x should be chosen so that it 
obeys Eq. 15. It is easy to see that the following change 
of variables does the job: 

^rn — ^ 2 m+l C2m—Ij ^ — I5 ■■■5 Ef 1, 

C'n = 2'KQj -|- 2ttQ\^, ^'n+1 ~ ^ 1 ) 

= C2m-2N-2 — C 2 N, m = N + 2 ,..., 2N 
^2N+1 ~ C 2 N — Cl — 2'KkQ^+TT^ C2^_|_2 = —27r(5^, 


Indeed, it is easy to check that 

_ /Oat+i —V 
\ E) Ojv+i 

where T) is the -|- 1 dimensional diagonal matrix 


(2k 0 ••• 0 0 \ 
0 2k ■■■ 0 0 


V = 


0 0 ••• 2 0 

Vo 0 ••• 0 1/ 


We can see that this is in the canonical skew-normal form 
shown in Eq. 14, with the parameters M = 2N + 2, 
7 = iV -h 1 and 


di = ... = djM—i = 2fc, cZat = 2, dAf-i-i = 1 

With these results we can write down the low energy 
effective Hamiltonian Hf^g: according to Eq. 16, H^g 
must take the form 

2N 

77eff = EE vpal^apm (65) 

m—l p>0 

where the sum runs over p = 7 r/ 2 s, 37 r/ 2 s,.... Notice that 
Heg does not include a term of the form ■.., C^j) 

which was present in the previous examples. The reason 
that this term is not present is that M = 27 in this case 
— that is, none of the C' terms commute with all the 
other Cj. This is closely related to the fact that the mo¬ 
mentum p is quantized in half-odd integer multiples of 
7r/2s so unlike the previous examples, we cannot con¬ 
struct an operator Upm with p = 0 (sometimes called a 
“zero mode” operator). 

Let us now discuss the effective Hilbert space TZeff- 
According to the results of section HB, the effective 
Hilbert space HeS is spanned by states {\a,{npm})} 
where \cy.,{npm}) is the unique simultaneous eigenstate 
of the form 

g*C''/ 2 fc|^,{^^^}) ^ - 1 

E^'«/'^\oL,{np^}) = e”“”|a,{np„}), 
e |q:, {ripm}) = |q:, {ripm}), 

E^'\cx,{npm}) = \cx,{npm}), i = N + 2, ...,2N + 2 

1^5 I {^pm}) (^ 6 ) 

Here the label Upm runs over non-negative integers, while 
ot is an abbreviation for the N component integer vector 
(ai, ...,aN) where runs over two values { 0 , 1 }, and 
the other a^’s run over ( 0 , 1 ..., 2k — 1 }. 

As in the previous cases, we can easily diagonalize the 
effective theory: clearly the {|q:, {ripm})} basis states are 
also eigenstates of HeS with energies given by 

N 

^=EE vpnpm (67) 

m—l p>0 
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The above equation gives the complete low energy spec¬ 
trum of H in the limit C/ —> oo. 

An important feature of the above energy spectrum 
(67) is that the energy E is independent of a. It fol¬ 
lows that every state, including the ground state, has a 
degeneracy of 

D = 2- (68) 

since this is the number of different values that a ranges 
over. 

We now discuss the physical meaning of this degener¬ 
acy. As in the previous examples, when [/ —>■ oo, the im¬ 
purities reconnect the edge modes, breaking the edge up 
into 2N disconnected components associated with the in¬ 
tervals [x 2 n,Xi],[xi,X 2 ],.:,[x 2 n-i,X 2 n] (Fig- 4b). The 
npm quantum numbers describe the number of phonon 
excitations of momentum p in the mth component of the 
edge. The a quantum numbers also have a simple physi¬ 
cal interpretation. Indeed, if we examine the definition of 
OL (66), we can see that for i ^ N, where 

Qi is the total charge in the interval [x 2 i-i,X 2 i+i\ while 
giTrajv _ gjTrq q jg ^j^g total charge on the edge. 

Thus, for i ^ N, ai/k is the total charge in the interval 
[x 2 i-i, X 2 i+i] modulo 2 while aN is the total charge on 
the edge modulo 2. The quantum number ranges over 
two possible values {0,1} since the total charge on the 
edge must be an integer while the other ai’s range over 
2k values {0,1,..., 2k — 1} since the fractional quantum 
spin Hall state supports excitations with charge 1/k and 
hence the charge in the interval [x 2 i-i, X 2 i+i] can be any 
integer multiple of this elementary value. 

It is interesting to compare our formula for the degen¬ 
eracy (68) to that of Refs. 23-27. Those papers studied a 
closely related system consisting of a FQSH edge in prox¬ 
imity to alternating ferromagnetic and superconducting 
strips (Fig. 5a). The authors found that this related 
system has a ground state degeneracy of 13 = {2k)^~^, 
which agrees exactly with our result, since Refs. 23-27 
did not include the two-fold degeneracy associated with 
fermion parity. In fact, it is not surprising that the two 
systems share the same degeneracy since one can tune 
from their system to our system by shrinking the size 
of the ferromagnetic and superconducting strips while at 
the same time increasing the strength of the proximity 
coupling (see Fig. 5b-c). 

Although our system shares the same degeneracy as 
the one studied in Refs. 23-27, one should keep in mind 
that there is an important difference between the two de¬ 
generacies: the degeneracy in Refs. 23-27 is topologically 
protected and cannot be split by any local perturbation, 
while our degeneracy is not protected and splits at any 
finite value of U, as we explain in the next section. That 
being said, if we modify our model in a simple way, we can 
capture the physics of a topologically protected degener¬ 
acy. In particular, the only modification we would need 
to make is to replace each individual magnetic impurity 
with a long array of many magnetic impurities, and sim¬ 
ilarly we would replace each individual superconducting 



FIG. 5. (a) A fractional quantum spin Hall edge in proxim¬ 
ity to alternating ferromagnetic and superconducting strips, 
(b-c) By shrinking the size of the ferromagnetic and supercon¬ 
ducting strips while at the same time increasing the strength 
of the proximity coupling, we can continuously deform the 
system into a fractional quantum spin Hall edge with mag¬ 
netic and superconducting impurities. 


impurity with a long array of many superconducting im¬ 
purities. After making this change, the degeneracy would 
remain nearly exact even at finite U, with a splitting 
which is exponentially small in the length of the arrays. 


E. Finite U corrections 

In the previous sections we analyzed the low energy 
physics of three different systems in the limit t/ —>■ oo. 
In this section, we discuss how these results change when 
U is large but finite. 


1. Single magnetic impurity 


We begin with the simplest case: a single magnetic im¬ 
purity on the edge oi & v = 1/k fractional quantum spin 
Hall state. We wish to understand how finite U correc¬ 
tions affect the low energy spectrum derived in section 
HIB. 

We follow the approach outlined in section HD. Ac¬ 
cording to this approach, the first step is to construct 
the operator H which is conjugate to the argument of the 
cosine term, C = fc((/i^(0)-|-<)>^(0)). To do this, we regular- 

ize C as in equation (40), replacing C —>■ /_^^2 + 

(j)^{x))5{x)dx. For concreteness, we choose the regulated 
delta function 5 to be 


5{x) 


i |a;| < b/2 
0 |a;| > b/2 


With this regularization, we find: 


= -y 


kvi 


[C, [C,Ho]] = - 


-b/2 

Airkv 


{dxi/t 


dj;(j)i)dx 
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so that 


ward calculation gives 


1 

n = ^ / - dx(l)i)dx 

^ J-b/2 

According to Eq. 23, the low energy theory at fi¬ 
nite U is obtained by adding terms to iJeff (46) of the 
form YlnZ-oc e“^e„({op, a],}, Cg). Here, the e„ are some 
unknown functions whose precise form cannot be deter¬ 
mined without more calculation. We should mention that 
the e„ functions also depend on ?7 — in fact, Cn —t 0 as 
17 —>■ oo — but for notational simplicity we have chosen 
not to show this dependence explicitly. In what follows, 
instead of computing e„, we take a more qualitative ap¬ 
proach: we simply assume that e„ contains all combina¬ 
tions of Up, aj), Cg that are not forbidden by locality or 
other general principles, and we derive the consequences 
of this assumption. 

The next step is to analyze the effect of the above terms 
on the low energy spectrum. This analysis depends on 
which parameter regime one wishes to consider; here, we 
focus on the limit where L —)■ oo, while U is fixed but 
large. In this case, HeS (46) has a gapless spectrum, so we 
cannot use conventional perturbation theory to analyze 
the effect of the finite U corrections; instead we need to 
use a renormalization group (RG) approach. This RG 
analysis has been carried out previously^^42 
not repeat it here. Instead, we merely summarize a few 
key results: first, one of the terms generated by finite U, 
namely is relevant for fc > 1 and marginal for k = 1. 
Second, the operator can be interpreted physically as 
describing a quasiparticle tunneling event where a charge 
1 /fc quasiparticle tunnels from one side of the impurity 
to the other. Third, this operator drives the system from 
the U = oo fixed point to the U = 0 fixed point. 

These results imply that when fc > 1, for any finite 
[7, the low energy spectrum in the thermodynamic limit 
L — )■ 00 is always described by the U = 0 theory Hq. 
Thus, in this case, the finite U corrections have an im¬ 
portant effect on the low energy physics. We note that 
these conclusions are consistent with the RG analysis of 
magnetic impurities given in Ref. 32. 


2. Multiple magnetic impurities 


We now move on to consider a system of N equally 
spaced magnetic impurities on an edge of circumference 
L. As in the single impurity case, the first step in un¬ 
derstanding the finite U corrections is to compute the H^ 
operators that are conjugate to the Q’s. Regularizing 
the cosine terms as in the previous case, a straightfor- 


Aiij 




nb 

1 
2 


i+b/2 


Xi-b/2 




dx(l)\.)dx 


where i,j = 1 ,...,A^.^^ 

According to Eq. 23, the finite U corrections 
contribute additional terms to H^g (54) of the form 

where n = (ni,...,nAr) is an N- 
component integer vector. Here, Cn is some unknown 
function of the operators {opm, a|,m) ^m} which vanishes 
as t/ —>■ oo. 

We now discuss how the addition of these terms af¬ 
fects the low energy spectrum in two different parameter 
regimes. First we consider the limit where L ^ oo with 
U and N fixed. This case is a simple generalization of the 
single impurity system discussed above, and it is easy to 
see that the same renormalization group analysis applies 
here. Thus, in this limit the finite U corrections have a 
dramatic effect for fc > 1 and cause the low energy spec¬ 
trum to revert back to the U — 0 system for any finite 
value of 17, no matter how large. 

The second parameter regime that we consider is where 
L,N —i' oo with U and L/N fixed. The case is different 
from the previous one because 77eff (54) has a finite en¬ 
ergy gap in this limit (of order v/s where s = L/N). 
Furthermore, iJeff has a unique ground state. These 
two properties are stable to small perturbations, so we 
conclude that the system will continue to have a unique 
ground state and an energy gap for finite but sufficiently 
large U. 

The presence of this energy gap at large U is not sur¬ 
prising. Indeed, in the above limit, our system can be 
thought of as a toy model for a fractional quantum spin 
Hall edge that is proximity coupled to a ferromagnetic 
strip. It is well known that a ferromagnet can open up 
an energy gap at the edge if the coupling is sufficiently 
strong^^’^^, which is exactly what we have found here. 


3. Multiple magnetic and superconducting impurities 


Finally, let us discuss a system of 2N equally spaced 
alternating magnetic and superconducting impurities on 
an edge of circumference L. As in the previous cases, the 
first step in understanding the finite U corrections is to 
compute the H^ operators that are conjugate to the Ci’s. 
Regularizing the cosine terms as in the previous cases, a 
straightforward calculation gives 




H, = 


1 

2 


f>Xi-\-bj2 
ri — bj2 


- {-^y^^dx(j)i)dx 


where i,j = 1,..., 2N. 
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As a first step towards understanding the finite U cor¬ 
rections, we consider a scenario in which only one of the 
impurities/cosine terms has a finite coupling constant t/, 
while the others have a coupling constant which is in¬ 
finitely large. This scenario is easy to analyze because we 
only have to include the corrections associated with a sin¬ 
gle impurity. For concreteness, we assume that the impu¬ 
rity in question is superconducting rather than magnetic 
and we label the corresponding cosine term by cos(C' 2 j) 
(in our notation the superconducting impurities are la¬ 
beled by even integers). Having made these choices, we 
can immediately write down the finite U corrections: ac¬ 
cording to Eq. 23, these corrections take the form 

OO 

n——oo 

where the e„ are some unknown functions which vanish 
as t/ — >■ OO. 

Our next task is to understand how these corrections 
affect the low energy spectrum. The answer to this ques¬ 
tion depends on which parameter regime one wishes to 
study: here we will focus on the regime where L, > oo 
with L/N and U fixed. In this limit, (65) has a 
finite energy gap of order v/s where s = L/N. At the 
same time, the ground state is highly degenerate: in fact, 
the degeneracy is exponentially large in the system size, 
growing as D = 2- (2/c)^“^. Given this energy spectrum, 
it follows that at the lowest energy scales, the only effect 
of the hnite U corrections is to split the ground state 
degeneracy. 

To analyze this splitting, we need to compute the ma¬ 
trix elements of the finite U corrections between different 
ground states and then diagonalize the resulting D x D 
matrix. Our strategy will be to use the identity (26) 
which relates the matrix elements of the finite U correc¬ 
tions to the matrix elements of . Following this 

approach, the first step in our calculation is to compute 
T 2 j. Using the definition (24), we find 

r 2 i = - ^2^-1) 

assuming j ^ N. (The case j = N is slightly more com¬ 
plicated due to our conventions for describing the peri¬ 
odic boundary conditions at the edge, so we will assume 
j 7 ^ A in what follows). 

The next step is to find the matrix elements of the 
operator between different ground states. To this 

end, we rewrite T 2 j in terms of the C' operators: T 2 j = 
c' 

. The matrix elements of can now be computed 

straightforwardly using the known matrix elements of C/ 
(see Eqs. D7-D8): 

la) = e ^ Sa'a. (69) 

where |a) denotes the ground state |a) = |a,n = 0 ). 

At this point we apply the identity (26) which states 
that the matrix elements of e®”^^^e„({apm, «pm}) 


equal to the matrix elements of Un ■ where is 

some unknown proportionality constant. Using this iden¬ 
tity together with (69), we conclude that the matrix ele¬ 
ments of the finite U corrections are given by 

where f{x) = J2'^=-oo Une™^. 

We are now in position to determine the splitting of the 
ground states. To do this, we note that while we don’t 
know the values of the constants and therefore we 
don’t know f{x), we expect that generically the function 
/ will have a unique minimum for aj S { 0 ,1,..., 2k — 1 }. 
Assuming this is the case, we conclude that the Hnite U 
corrections favor a particular value of say aj = 0 . 
Thus these corrections reduce the ground state degener¬ 
acy from D to D/2k. 

So far we have analyzed the case where one of the su¬ 
perconducting impurities is characterized by a finite cou¬ 
pling constant U while the other impurities are at infinite 
coupling. Next, suppose that all the superconducting im¬ 
purities have finite U while all the magnetic impurities 
have infinite U. In this case, similar analysis as above 
shows that the matrix elements of the finite U correc¬ 
tions are of the following form:^^ 



To determine the splitting of the ground states, we need 
to understand the eigenvalue spectrum of the above ma¬ 
trix. Let us assume that / has a unique minimum at 
some aj = q — which is what we expect generically. 
Then, as long as the system size is commensurate with 
this value in the sense that Nq is a multiple of k, we 
can see that the above matrix has a unique ground state 
|a) with tti = ... = oat-i = q and a^ = Nq/k (mod 
2). Eurthermore, this ground state is separated from the 
lowest excited states by a finite gap, which is set by the 
function /. Thus, in this case, the finite U corrections 
completely split the ground state degeneracy leading to 
a unique ground state with an energy gap. 

Likewise, we can consider the opposite scenario where 
the magnetic impurities have finite U while the supercon¬ 
ducting impurities have infinite U. Again, similar anal¬ 
ysis shows that the corrections favor a unique ground 
state which is separated from the excited states by a fi¬ 
nite gap. The main difference from the previous case is 
that the matrix elements of the finite U corrections are 
off-diagonal in the |a) basis so the ground state is a su¬ 
perposition of many different |a) states. 

To complete the discussion, let us consider the case 
where all the impurities, magnetic and superconducting, 
have finite U. If the magnetic impurities are at much 
stronger coupling than the superconducting impurities or 
vice-versa then presumably the finite U corrections drive 
the system to one of the two gapped phases discussed 
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above. On the other hand, if the two types of impurities 
have comparable values of U, then the low energy physics 
is more delicate since the finite U corrections associated 
with the two types of impurities do not commute with 
one other, i.e. ^ 0- In this case, a more 

quantitative analysis is required to determine the fate of 
the low energy spectrum. 


IV. CONCLUSION 

In this paper we have presented a general recipe for 
computing the low energy spectrum of Hamiltonians of 
the form (1) in the limit U ^ oo. This recipe is based 
on the construction of an effective Hamiltonian H^g and 
an effective Hilbert space HeS describing the low energy 
properties of our system in the infinite U limit. The 
key reason that our approach works is that this effective 
Hamiltonian is quadratic, so there is a simple procedure 
for diagonalizing it. 

While our recipe gives exact results in the infinite U 
limit, it provides only approximate results when U is fi¬ 
nite; in order to obtain the exact spectrum in the finite U 
case, we need to include additional (non-quadratic) terms 
in Hf,g. As part of this work, we have discussed the gen¬ 
eral form of these finite U corrections and how they scale 
with U. However, we have not discussed how to actually 
compute these corrections. One direction for future re¬ 
search would be to develop quantitative approaches for 
obtaining these corrections — for example using the in- 
stanton approach outlined in Ref. 35. 

Some of the most promising directions for future work 
involve applications of our formalism to different physical 
systems. In this paper, we have focused on the applica¬ 
tion to Abelian fractional quantum Hall edges, but there 
are several other systems where our formalism could be 
useful. For example, it would be interesting to apply 
our methods to superconducting circuits — quantum cir¬ 
cuits built out of inductors, capacitors, and Josephson 
junctions. In particular, several authors have identified 
superconducting circuits with protected ground state de¬ 
generacies that could be used as qubits.The for¬ 
malism developed here might be useful for hnding other 
circuits with protected degeneracies. 
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Appendix A: Derivation of low energy effective 
theory 

In this appendix, we derive an effective theory that 
describes the low energy spectrum of Lf (1) in the limit 
U —>■ oo. More specifically, we show that the low energy 
spectrum of H in the infinite U limit is described by the 
effective Hamiltonian H^g (2), which is defined within the 
effective Hilbert space Heg (5). Before proving this result 
in generality, we first derive it for two illustrative exam¬ 
ples in appendix A 1 and a special case in appendix A 2. 
Finally, after this preparation, we work out the general 
case in appendix A 3. 

1. Two examples 

a. Harmonic oscillator with a cosine term 

To understand the basic ideas underlying the deriva¬ 
tion, it is helpful to consider some simple examples. We 
start by studying a one dimensional harmonic oscillator 
with a cosine term; 

Kx^ 

H = ^ ^ -U cos(27rx) (Al) 

2m 2 

In the following, we derive an effective Hamiltonian H^g 
and effective Hilbert space T-Leg that describe the low en¬ 
ergy spectrum of H in the limit U ^ oo. 

To begin, we decompose H into two pieces, H = Hi + 
H 2 , where 

Hi = - U cos(27rx), H 2 = —— (A2) 

2 m 2 

Our strategy is as follows: first, we show that when U 
is large. Hi has a collection of nearly degenerate ground 
states which are separated from the lowest excited states 
by a large gap. Next we argue that we can treat H 2 as 
a perturbation which splits the ground state degeneracy 
of Hi. Finally, using degenerate perturbation theory, we 
derive a low energy effective Hamiltonian for our system. 

Following this plan, we start with the Hamiltonian Hi. 
This Hamiltonian describes a one dimensional particle 
moving in a cosine potential. The low energy physics of 
Hi is especially simple when U is large. In this case, tun¬ 
neling between the different cosine minima is suppressed 
so that Hi has an infinite set of nearly degenerate ground 
states — one for each cosine minimum. We will label 
these states as {\ipq)} where is localized around the 
minimum x = q and <7 = 0, ±1, ±2,.... 

We can estimate the energy gap of Hi by expanding the 
cosine potential to quadratic order in x: U cos(27rx) « 
t7(l —27r^x^). In this approximation, the cosine potential 
is equivalent to a harmonic oscillator with frequency w = 
2'it^JU jm. In particular, it follows that the ground states 
of Hi are separated from the lowest excited states by an 
energy gap A of order A ~ i/ujm. 
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Now let us imagine adding H 2 to Hi. We would like to 
know how H 2 splits the degeneracy between the ground 
states \ipq). To answer this question, we will treat H 2 as 
a perturbation to Hi and then we will compute the as¬ 
sociated energy splitting using degenerate perturbation 
theory. However, before we do this, we first need to 
check that such a perturbative approach is justified in 
the infinite U limit. To this end, we need to estimate 
the size of the matrix element {ipex\H 2 \ipq) where j'i/'g) is 
an arbitrary ground state and \ipex) is an arbitrary ex¬ 
cited state of Hi. For large U we can approximate \tpq) 
by a harmonic oscillator ground state centered at posi¬ 
tion X = q. Similarly, we can approximate \ipex) by the 
nth excited state of the harmonic oscillator centered at 
X = q. Within this approximation, the matrix element 
{tpex\H 2 \'tpq) reduces to 

{'>Pex\H 2 \'ipq) ~ {n\K{x + g)^|0) (A3) 

where | 0 ) and |n) are the ground state and nth ex¬ 
cited state of a harmonic oscillator with frequency w = 
2'K^JTJJ^n, centered at a; = 0. The latter matrix element 
can be evaluated easily with the result 

f {muy/* if = 1 

{n\K{x + q)'^\0) ^ if n = 2 

/ I / I y/YYlU 

I 0 if n > 3 


Combining this expression with our formula for the en¬ 
ergy gap A, we obtain 


E 

Ibex) 


\{'ljjex\H2\'lpq)\^ 

A 


u 


(A4) 


in the large U limit. This estimate is significant because 
the left hand side of (A4) is proportional to the second 
order perturbative correction to the ground state ener¬ 
gies. Evidently, this correction vanishes as U —>■ 00 , so we 
conclude that first order perturbation theory gives exact 
results in this limit. 

With this justification, we now proceed with the per¬ 
turbative calculation. According to first order degenerate 
perturbation theory, the energy splitting of the ground 
states can be determined by diagonalizing the matrix 
{tpqi\H 2 \'ipq). When U is large, Itpq) can be approximated 
as a Gaussian wave function, centered at x = q. The 
width of this Gaussian is given by: 


(Ax)^ 


1 

muj 



We see that as U —>■ 00 , (Ax)^ —0 so that \ipq) ap¬ 
proaches a position eigenstate: IV'g) —5^ \x = q). We con¬ 
clude that the low energy spectrum of H can be obtained 
by diagonalizing the matrix {q'\H 2 \q). 

At this point, our calculation is essentially complete: 
the matrix elements {q'\H 2 \q) define our low energy ef¬ 
fective Hamiltonian, while the ground state subspace 


spanned by Hg)} defines our low energy Hilbert space. 
In other words, the low energy effective Hamiltonian is 
given by 

Kx^ 

H,s = H2 = ^ (A5) 

while the low energy effective Hilbert space HeS is the 
subspace spanned by position eigenstates Hg)} for q = 
0, ±1,±2.... Glearly this effective theory is valid for en¬ 
ergies smaller than the gap of Hi, i.e. E < 2-KyJU/m. 

Gomparing these expressions with the effective Hamil¬ 
tonian and Hilbert space from section HA, we see that 
they agree exactly. Thus we have successfully established 
equations (2) and (5) for this example. 


b. Harmonic oscillator with 2 cosine terms 


Another important illustrative example is given by 
a one dimensional harmonic oscillator with two cosine 
terms: 

Kx'^ 

H = -1- U cos(dp) — U cos(27ra:) (A6) 

2m 2 

Here, d is a positive integer. 

Before analyzing H we need to choose an appropriate 
basis in which to represent it. Because the arguments of 
the two cosine terms don’t commute with one another, 
neither the position basis nor the momentum basis are 
particularly convenient choices. Instead, we find it help¬ 
ful to work in a third basis, which consists of simulta¬ 
neous eigenstates of the commuting operators 
We will denote these simultaneous eigenstates by \9,ip) 
where 


E^\e,p) = E%e,p) 

e^-^^\e,p) = E^\e,p) (A7) 


Here, the labels 9,tp take values in 0 < 0,^5 < 2.tt. The 
explicit formula for \9,Lp) is 


\Q,'y) = 



(AS) 


where \ j + ^) denotes the position eigenstate at position 

We now work out what the Hamiltonian H looks like 
in the 9, cp representation. The first step is to express the 
x,p operators in terms of 9,ip. To this end, we observe 
that 

e^‘^^\9,ip) = e“‘^A^|6»-ha,v5) 
e“^|d, (fi) = \9,tf — 27ra) 


Differentiating these equations with respect to a, we de¬ 
rive 


x\9,p) 

p\9,p) 


1_9 ^ 

i 89 2 tt 


\S,^) 


27r 8 
i 8ip 
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From these equations, we deduce that 

{e,ip\p\ii) = 

for any state |'0)- We conclude that in the 0,Lp represen¬ 
tation, the operators x,p take the form 

19 If 2 tt d 

or equivalently 

X = -pe + P = ‘^'xp^ (A9) 

ZTT 

where pg = and p^ = 

The next step is to find expressions for 6^^^“ and 
in terms of 0, f. One way to do this is to exponentiate 
(A9): 

^2Trix _ ^-2TTipe+i(p ^ip _ ^2Trip^ 

We then make use of two operator identities which we 
will prove shortly: 

= e®® (AlO) 

With these identities, we can simplify the expressions for 

^2mx gip 

6 ^®®“ = e®"^, e®P = e®® (All) 

(Alternatively, we could also have derived (All) directly 
from (A7)). The proof of the identities (AlO) relies on 
the following observations: 

\e, f) = \e- 27r, f), |d, f) = \6,f- 27r) 

and 

\0 - 2tt, f) = \0, f), \0, f-2TT)= e®®|6», f) 


the boundary conditions for these wave functions, we use 
the identities 

|0,2^)=e-®«|0,O), \2Ti,f) = \Q,f) 

which follow from the definition of \0,f). These identities 
imply that our wave functions i/' satisfy the boundary 
conditions 

V’(6',27r) = e®®?/'(^',0), 'ip{2Tr,f) = 'ip{0,f) (A13) 


So far, all we have done is derive the 0, f represen¬ 
tation (A12) of the Hamiltonian and the Hilbert space 
(A13). We now use this representation to find the low 
energy spectrum of H in the large U limit. To begin, we 
note that H can be thought of as describing a particle on 
a torus parameterized by 0, f. This particle is coupled to 
a vector potential and two cosine potentials. Now con¬ 
sider the limit where U is large. In this limit, tunneling 
between the different cosine minima is suppressed, so we 
conclude that H has a set of nearly degenerate ground 
states, each of which is localized in a different minimum. 
There are d different cosine minima located at positions 
{0, f) = (^-Kajd, 0), with a = 0,1,..., d — 1, so H has d 
ground states. We will label these states by \4>a)- 

To estimate the energy gap separating the ground 
states from the lowest excited states, we expand the co¬ 
sine potentials to quadratic order. In this approximation, 
H reduces to a sum of two decoupled harmonic oscilla¬ 
tors with frequencies y^U/m and VUK. We conclude 
that the energy gap is of order A ^ min(-y/CZ/m, VUK). 

Let us now translate these results into the language 
of effective Hamiltonians. We have seen that H has d 
ground states \tjja)- We have also seen that these states 
are separated from the excited states by a large gap A. 
Furthermore, it is easy to see that when 17 —> oo, the 
width of |^/>a) becomes vanishingly small, so that {ipa) 
approaches the state \0 = ^^,f = 0): 


Putting these together, we deduce that 


g2«p«|0,^) ^ e^-^P-\0,f)=e^<>\0,f) 


\4>a) —t \0 


2'Ka 


f = 0) 


Since these relations hold for all basis states, they imply 
the operator identities (AlO). 

We now have all the ingredients to write the Hamilto¬ 
nian H in terms of 0,f'. combining (A9) and (All), we 
derive 

^ ^ 47 r^Py K{pg - fl2TTf 
2m 2 

— U cos[d0) — U cos{f) (A12) 

This Hamiltonian is defined on a Hilbert space consisting 
of wave functions '!/'(d, f) with Q < d^f < 27r. To find 


Putting this all together, we conclude that the low energy 
spectrum of H is described by an effective Hamiltonian, 
i7eff = 0 defined within an effective d-dimensional Hilbert 
space spanned by the states \0 = ^^,f = 0), with 
a = 0,1,d — 1. This effective description is valid for 
energies E < A. (Here, the reason we use i7eff = 0 rather 
than iJeff = const, is that we only interested in energy dif¬ 
ferences and therefore we are free to redefine the ground 
state energy to be 0). Comparing these results with the 
effective Hamiltonian and Hilbert space from section H A, 
we can see that there is exact agreement. 
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2. Special case 

We now generalize the example of appendix A 1 a to a 
Hamiltonian of the form 

M 

H = Hq — U'^^cos{2Trxi), (A14) 

i-l 

defined on the 2N dimensional phase space 
{xi,...,xn,Pi,-,Pn} with [xi,pj] = iSij. Here Hq 
is an arbitrary positive semidefinite quadratic function 
of {xi, ...,Xn,Pi, ■■■,Pn} with the only restriction being 
that the M x M matrix 

Afij = -[xi, [xj,Ho]] (A15) 


is non-degenerate. 

Following the same outline as in the previous sections, 
we first derive an effective Hamiltonian and effective 
Hilbert space that describe the low energy spectrum of 
H in the infinite U limit, and then we show that this 
effective theory agrees with the general expressions from 
equations (2) and (5). 

For the first step, we use the same strategy as in ap¬ 
pendix Ala: we decompose the Hamiltonian into two 
pieces, H = Hi + H 2 , where 


M M 

= E o '' • - t/^ cos(2^a:,) 


i,j = l 


M 




i,i=l 


2 ■ niHj 


(A16) 


Here A4ij is an M x M scalar matrix defined by A1 = 
J\f~^ and Hi, ...jHm are operators defined by 

M 

Ui =-i'^Mij[xj,Ho] (A17) 

1=1 

After making this decomposition, we will treat H 2 as a 
perturbation to Hi and then derive an effective Hamilto¬ 
nian using first order degenerate perturbation theory. 

Before executing this plan, we first make some prelim¬ 
inary observations. One observation is that 

[x^,Iij\= i5ij, z = l,...,M (A18) 

Another observation is that we can assume without loss 
of generality that 

[xj.Hj] = [pi.Hj] = 0, i = M +l,...,N (A19) 

The reason why we can assume (A19) is that we can 
always redefine the position and momentum operators 
Xi,Pi for j = M -I- 1, ...,N according to 

M M 

Xj -)■ Xj + ^ a^kXk, Pi ^ Pz + ^ bikXk 


where Uik = i[xi,Iik] and bik = i[pi^Iik]- After this re¬ 
definition, Eq. A19 is automatically satisfied. A final 
observation is that H^ can be written in the form 

= Pj + (A20) 

where Aj is a linear combination of {xi, ..., xm}- Indeed, 
this result follows immediately from (AI 8 ) and (A19). 

With these observations in mind, we now study the 
low energy spectrum of Hi. To begin we note that (A20) 
implies that Hi can be written in the form 


Hi 


M 


E 


2 


M 

■ {pi + Ai){pj +Aj)-U'Y^ cos(27rxi) 
i=l 


where Aj is a linear function of {xi, ...,xm}- Next, we 
note that [xi, Hi] = 0 for i = M + 1 ,..., N, which implies 
that Hi can be diagonalized separately for each value 
of x±_ = {xm+1, Once we fix x±, the Hamilto¬ 

nian Hi describes an M dimensional particle with co¬ 
ordinates {xi, ...,xm), moving in a periodic potential 

— U X]i=i cos(a;i) and coupled to a vector potential A that 
depends linearly on {xi, ..., xm}- Let us consider the low 
energy physics of this M dimensional particle when U is 
large. In this case, we can neglect tunneling between the 
different minima of the cosine potential, and treat each 
minimum in isolation. At the same time, it is easy to 
see that the energy spectra associated with different co¬ 
sine minima and different values of a;_L are all identical 
since Hi has discrete (magnetic) translational invariance 
in the xi,..., xm directions as well as continuous transla¬ 
tional invariance in the xm+i, ■■■,Xn directions. Putting 
these facts together, we conclude that the Hamiltonian 
Hi has an infinite set of nearly degenerate ground states 

— one ground state for each cosine minimum and each 
value of x±. We label these ground states as \tpq.x±) 
where q = {qi, ...,qM) is an M component integer vec¬ 
tor describing the position of the cosine minimum and 
xx = [xm+i, ■■■,xn) is an N — M component real vector 
describing the position in the orthogonal directions. 

We can estimate the energy gap between the ground 
states and excited states of Hi by expanding the co¬ 
sine potential to quadratic order in a:. In this approx¬ 
imation, the cosine potential reduces to a multidimen¬ 
sional quadratic potential. Diagonalizing this potential 
gives a collect ion of harmonic oscillators with frequencies 
LOi = \JWJmi where are the eigenvalues of the matrix 
M-ij. The energy gap is determined by the smallest fre¬ 
quency and thus the largest eigenvalue We conclude 
that Hi has an energy gap A ~ jm where m is the 
maximum eigenvalue of Mij. 

We now imagine adding H 2 to iJi, and we ask how 
the low energy spectrum changes. As in appendix Ala, 
we will answer this question using a perturbative expan¬ 
sion in H 2 . However, before we perform this calcula¬ 
tion, we need to check that this perturbative approach 
is valid when U is large. To this end, we need to es¬ 
timate the size of the matrix elements {'il’ex\H 2 \ilJq,x±) 
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where IV^q^cc^), \'4>ex) are arbitrary ground states and ex¬ 
cited states of Hi. The first step is to observe that H 2 
has a special property: for any z = 1 ,M, we have 

M 

[xi, H2] = [xi, Hq] - i ■ n^- 

= [xi,HQ] - [x„Hq] 

= 0 

It follows that the momentum operators pi,...,pM do 
not appear in 7 ? 2 - Thus, H2 is a sum of three types of 
terms: XiXj, pkPi and XiPu, where i,j are arbitrary and 
k,l > M + 1. We need to estimate the matrix elements 
corresponding to each of these terms. We can do this 
using the same argument as in appendix Ala. First, we 
note that when U is large, the states \'ipq.xx_) and IV'ea;) 
can be approximated as harmonic oscillator eigenstates 
centered at some appropriate positions in space. The ma¬ 
trix elements of interest can then be related to harmonic 
oscillator matrix elements as in Eq. A 3 . Omitting de¬ 
tails, a straightforward calculation shows that all three 
types of matrix elements fall off like or faster as 

t/ —^ 00. Combining this scaling law with the expression 
for A, we see that 

\{4’ex\H2\tpq)\'^ ^ Q f 

IV’ex) 

As in appendix Ala, this estimate implies that the sec¬ 
ond order perturbative corrections to the ground state 
energies vanish in the limit 1 / —> 00. Therefore, first 
order perturbation theory is exact in this limit. 

We now proceed with the perturbative calculation. Ac¬ 
cording to first order degenerate perturbation theory, the 
energy splitting of the ground states can be determined 
by diagonalizing the matrix {'4’q',x'^\H2\i^q.x±)- These 
matrix elements are easy to compute. Indeed, for large 
U, the ground states IV'q.cc^) can be approximated as 
harmonic oscillator ground states with a width of or¬ 
der ~ A-ft - In the limit U —)■ 00, the width 

' ' y/mU 

Ax —0 so |' 0 q,xj_) approaches a position eigenstate: 
IV'q.xx) |q, a:_L) where |q, a:_L) denotes the position 
eigenstate located at a: = {q,x±). Thus, in this limit, 
the matrix we need to diagonalize is {q',x'j_\H2\q,x±). 

Our calculation is now complete: we have shown that 
the low energy spectrum of H in the limit ?7 —>■ 00 can 
be obtained by diagonalizing the operator H 2 within the 
subspace spanned by {|q,a;j_)}. In other words, the ef¬ 
fective Hamiltonian for our system is 

HeS = H2 

M , .._i. 

= Ho-J2 2 '' ■ 

ij = l 

while the effective Hilbert space HeS is spanned by posi¬ 
tion eigenstates {|<7,a;j_)} for which the first M compo¬ 
nents qi,..., Qm are integers and the last N — M compo¬ 
nents xm+ 1 , ■■■,X]s[ are real valued. This effective theory 
is valid for energies E < 


Comparing the above effective Hamiltonian and 
Hilbert space with the general expressions from equa¬ 
tions (2) and (5), we can see that they match exactly. 
This completes our proof of equations (2), (5) for the 
special case (A14). 

3. General case 

To complete the derivation of equations (2) and (5), 
we now consider the most general case: 

M N 

H = Ho - cos(C'j), Ci = '^AijXj + AijPj + 

i=i j=i 

(A22) 

Here H is defined on a 2N dimensional phase space 
{xi,...,XnjPij--;Pn} with [x^.pj] = iSij. The first 
term Hq is a positive semidefinite quadratic function 
of {xi, ...,xn,Pi, ■■■,Pn} and 'jij.AijAi are arbitrary 
real numbers with the only constraints being that (1) 
{Cl, ...,Cm} are linearly independent, (2) [Ci,Cj] is an 
integer multiple of 27rz for all i,j so that the different 
cosine terms commute with one another, and (3) the ma¬ 
trix 

^v = -^[AACj.Ho]] (A23) 

is non-degenerate. Our task is to derive an effective 
Hamiltonian and effective Hilbert space that describes 
the low energy spectrum of H in the limit C/ —>■ 00. We 
will then show that this effective theory is exactly the 
one defined in equations (2) and (5). 

Our basic strategy is simple: we will map H onto the 
Hamiltonian studied in the previous section and then we 
will derive the effective theory using our previous results. 
While this is a simple plan at a conceptual level, there are 
some technical obstacles that make it difficult to define 
the desired mapping in one step. Therefore, we will build 
up the mapping in stages by making several successive 
changes of variables. 

In the first change of variables, we replace Xi,pi by new 
coordinates x[,p[ which are chosen so that the constraints 
Ci can be written as integer linear combinations of 27ra;' 
and p'. To this end, it is helpful to consider the M x M 
matrix Zij defined by 

Z.I = ^^[C.,C,] (A24) 

Clearly the matrix Zij is integer and skew-symmetric. 
Therefore, there exists an integer matrix V with deter¬ 
minant ±1 such that VZV^ = Z' where Z' is in skew- 
normal form^®: 
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Here 0 < I < M/2 and 0/ denotes an / x / matrix of 
zeros, and similarly for Om- 27- Let 

M 

= (A26) 

i=i 


with 

e*P^|0,V3,a;') = e*®’|0, yj, a;'), i = 1,...,/ 
i;'|0,y),a;') = x'/\9,ip,x'), i = / + 1,..., 


where 


Xi = TT • ^ Vy VjfeZjfc (mod 27r) (A27) 

j<k 

By construction, [C'',C''] = 2niZlj. Furthermore, be¬ 
cause we chose the offset Xi of fLe above form, we have 
the identity 


Here 0 < 6i, (fi < 27r while x'^ can be arbitrary real num¬ 
bers. The formal definition of \9,(p,x') is 


\9,ip,x') 


^ kiOi 

(2^)//2 ^ 

\ki + ..., k[ + ..., x'^^fASO) 


M 

(A28) 

i=i 

as one can check using the Campbell-Baker-Hausdorff 
formula. Equivalently, using the fact that Vij is an inte¬ 
ger matrix with determinant ±1, the above identity can 
be inverted and rewritten as 


where the ket on the right hand side denotes a position 
eigenstate localized at position = fci etc. 

We now re-express the Hamiltonian H in the 9, ip, x' 
representation. The first step is to find expressions for 
the x',p' operators in terms of 0i,ipi. Following the same 
logic as in appendix Alb, it is easy to show that x'^,p[ 
take the form 


M 

^iC, ^ giV-j^cr 
1=1 

The two identities (A28-A29) will be useful below. 

We are now in a position to construct the new variables 
x[,p[ that we seek. Specifically, we define 


, 19 <pi , 27r d 

i dOi 27r’ i d(pi 

for 7 = 1,...,/. Equivalently, we can write this as 


(A31) 


^'i = -Pei+’^^ Pz=‘^'^Pvi (A32) 


p-= for 1 = 1,...,/ 

= for7 = l,...,M-/ 

and we define p/^i ,..., p(v ^'m-i+i > • ■ • > ^'n f o be some 

arbitrary linear combination of Xi,pi with the only con¬ 
straint being that they obey the correct commutation 
relations. In the new variables, the Hamiltonian takes 
the form 

M / M 

// = //o - c/^cos - Xj] 

z=l \j=l 


Note that the argument of the cosines are now integer 
linear combinations of 27ra:' and p', just as we wanted. 

In the next step we find an alternative representation 
of H in which the arguments of the cosine terms all com¬ 
mute with each other. We accomplish this by working in 
an unusual basis, which is similar to the one discussed 
in appendix Alb. Specifically, let us consider the ba¬ 
sis of simultaneous eigenstates of the A^ -|- / commuting 
operators 


r iW-, ipi 27Tix'-, 27rixt f t i 

{e^i,...,e^Ee E...,e E Xj+i,..., a;^| 


We denote these simultaneous eigenstates by 


\9,ip,x') = \ei,..,ei,pi,...,ipi,x/^^,...,x'j^) 


where pe- = and p^p^ = Likewise, when i > 

/-I- 1, it is easy to see that a;',p' take the usual form, i.e. 
P^ = lT^p etc. 

In addition to x',p', we also need expressions for 
and for 1 = I,...,/. We can obtain these expressions 
by exponentiating (A32): 

g27rix' _ ^-2Tripg.+iipi _ ^2TTip^. 

As in appendix Alb, these equations can be simplified 
to 


^2TTix[ _ ^iipi ^ip[ _ giSi 

using the operator identities 

g27rip8. _ 2^ _ ^iOi 


(A33) 


(For an explanation of where these identities come from, 
see the discussion in appendix Alb). 

With this preparation, we are now ready to express the 
Hamiltonian H in terms of 9, tp, x'. Using (A33) together 
with the identity (A29), we can rewrite the cosine terms 
as 


M 


M 


=cos 27r^(V 


u=i 


i=i 


COS 
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where 


= 

/ di0i di9i 
I 27r ’■■■’ 27r ’ 27r’"'’ 


‘fl 

27r’ 


C/+1, 2^M-/ 


(A34) 


Similarly, using (A32) we can write the quadratic term 
Hq as a function of {9i, 6*/, x^} and 

{Pe^, Thus, all together, 

the Hamiltonian is given by 


M 

H = Hq — U ^ cos 

i=l 


M \ 

27tY,{V-%^A 


(A35) 


This Hamiltonian is defined on a Hilbert space consist¬ 
ing of wave functions ilj{9,ip,x') with 0 < 9i,ipi < 27r 
and arbitrary real numbers. As in appendix 

Alb, these wave functions obey boundary conditions of 
the form 


'fp\ei= 2 ^ = fp\ei=o for 1 = 1,...,/ (A36) 


and we define Xm+Ii ■■■Xn+i and pi, arbitrarily 

as long as they obey the canonical commutation relations. 
After this change of variables, our Hamiltonian takes the 
form 


M 

H ^ Ho-U'^cos (27rii) (A38) 

i=l 

At this point we have achieved the desired mapping: 
we can see that the above Hamiltonian (A38) is of the 
same form as (A14). This is very convenient because it 
means we can write down a low energy effective theory 
in the limit U ^ oo using our previous results without 
doing additional work. Indeed, according to Eq. A21, we 
have 


M / 1 _2 \ 

Heff = Ho-Yl 2 '' ■ 

where Aiij and H^ are defined by 


To proceed further we make an approximation in which 
we temporarily ignore the periodic boundary conditions 
on (pii9i. That is, we treat the Hamiltonian (A35) as if it 
were defined in a Hilbert space consisting of wave func¬ 
tions ip{6,(p,x') where 9i,p>i range over —oo < 9i,ipi < 
oo. Then, at the very end of our calculation we will rein¬ 
corporate the fact that 9i,ipi are actually angular vari¬ 
ables which range from 0 to 27r. The justification for this 
approximation is that the low energy eigenstates of H 
are localized near the minima of the cosine potential and 
are therefore insensitive to the global boundary condi¬ 
tions in the large U limit. Thus, the only effect of the 
periodic boundary conditions on the low energy theory is 
to identify certain states in the low energy Hilbert space 
■Heff — an effect we can take account of at the end of our 
derivation. 

Once we make this approximation, the variables 9^, ipi 
and x'j^-j^, ...jX'pf are all on an equal footing. Our 
Hamiltonian is then defined in a 2N + 21 dimen¬ 
sional phase space consisting of the (real valued) posi¬ 
tion variables {9i, ...,9j, (pi, ...ipj,x'j_^_^, ...,x'j^} together 
with their canonically conjugate (real valued) momenta 
{pe^,...,pej,p^i,...,p^i,p'i+i,...,p'N}- This completes 
the second step of our derivation: we have successfully 
rewritten the Hamiltonian in a form (A35) in which the 
arguments of the cosine terms commute with one an¬ 
other. 

In the final step, we make yet another change of vari¬ 
ables, defining new position operators ii,..., xat+j and 
new momenta Pi, ...,pn+i which are linear combinations 
of the previous position and momenta operators. The 
goal of this transformation is to simplify the arguments 
of the cosine terms even further. More specifically, we 
define 

M 


M=N Nij = -[xi, [xj,HQ]] 

and 

M 

H^ — i ^ {xj 7 Hq] 

i=i 

To complete the derivation, we need to express the 
effective Hamiltonian (A39) in terms of the original vari¬ 
ables Xi,pi. We will do this by relating the commutator 
[ii, Hq] to [Ci, Hq], To this end, we make a few observa¬ 
tions. First, we note that 

[2Trp^,-9„Ho] = 0 (A40) 

[pe„Ho] = 0 (A41) 

Here the first commutator vanishes because pi only ap¬ 

pears in Hq in the combination (—-|- ^Pi), while the 
second commutator vanishes because 9i doesn’t appear 
in Hq at all (see (A32)). Next, we observe that 

[9, ,Ho] = 27r [p^., Ho] = [p', Ho] (A42) 

where the first equality follows from (A40) and the second 
equality follows from (A32). Similarly, we have 

^[p^,Ho] = [-pe- + ^pi,Ho] = [x',Ho] (A43) 

where the first equality follows from (A41) and the second 
equality follows from (A32). If we now multiply (A43) 
by 27r and (A42) by di, we derive 

2^[C„Ho] = [C',Ho] 

for 1 = 1,..., M. It follows that 


i = (A37) 


27r[i„Ho] = [a. Ho] 


(A44) 
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With the help of (A44), it is now straightforward to 
rewrite i/eff in terms of Xi,pi. First, we take the commu¬ 
tator of (A44) with ij, and apply the Jacobi identity to 
derive AT = JV where Af is defined in Eq. A23. Likewise, 
we have A4 = A4 where A4 = Af~^. Substituting these 
relations into (A39), we see that our effective Hamilto¬ 
nian can be equivalently written as 


Appendix B: Diagonalizing the effective theory 

In this appendix, we derive a general recipe for diag¬ 
onalizing the effective Hamiltonian Heff (2). This recipe 
is the one outlined in section HB. Our analysis can be 
divided into three steps. 


iLeff = Ho-Y^ % . n,n 






where 


M 


“ 2Tri ^ 




1. Step 1: Creation and annihilation operators 

(A45) 

In general, the key to diagonalizing quadratic Hamil¬ 
tonians is to find appropriate creation and annihilation 
operators. For the case of JleB, this can be accomplished 
by searching for all operators a that are linear combina¬ 
tions of {xi, ...,X]v,pi, ...,p]v} and that satisfy 

[a, Heff] = Ea, (HI) 


Here we think of HojIIi, etc, as functions of the original 
variables Xi,pi. 

So far we have focused on the low energy effective 
Hamiltonian iLeffj but we also need to discuss the low 
energy effective Hilbert space 'Heff in which this Hamil¬ 
tonian is defined. If we apply the results of the previ¬ 
ous section, we see that T-LeS consists of all states sat¬ 
isfying Xi = (integer) for i = 1,...,M. Let us try to 
translate this back into our original variables. First, 
from Eq. A37 we see that HeS can be equivalently de¬ 
scribed as consisting of all states with = (integer) for 
i = 1,...,M (here we use the fact that Vij is an inte¬ 
ger matrix with determinant ±1). Next, we convert to 
the 9,(p,x' description, using the definition (A34). In 
the 6, ip, x' language, T-LeS consists of all states \6, p, x') 
with 9i = 2TTjdi • (integer), ipi = 27r ■ (integer), and 
x' = (integer) for i = / -|- 1,..., M — I. 

It is at this point we should remember that 9i and pi 
are actually angular variables (see the discussion below 
(A36)); that is we should identify 9i -|- 27r with 9i and 
Pi + 2 tt with Pi. Thus, the basis states \6,p,x') for 
the Hilbert space HeS can be parameterized in a non- 
redundant fashion by fixing pi = 0 and letting 9i range 
over the values 9i = 2'Kai/di where ai = 0,1, ...,di — 1. 
Meanwhile, since x' is not angular valued, it can range 
over arbitrary integers for i = I + 1,M — I and arbi¬ 
trary real numbers for j = M — / -|- 1,..., N. 

The above description in terms of the 9, p, x' variables 
is perhaps the most explicit way to parameterize the low 
energy effective Hilbert space 'HeS- However, it is also 
useful to describe "Heff in terms of our original variables 
Xi,pi. One way to do this is to note that the above 
basis states \9^p,x'_i_) are precisely the states that sat¬ 
isfy the constraints cos(C'i) = 1 for all i. Thus we can 
equivalently describe HeS as the set of states satisfying 
cos(C'i) = 1 for all i. 

This completes our derivation: we can see that our 
effective Hamiltonian iLefi and effective Hilbert space "Heff 
exactly match the expressions in equations (2) and (5). 
Thus we have derived these results in complete generality. 


for some scalar E ^ 0, as well as 

[a,Q]=0, z=l,...,M (B2) 

While the first condition (Bl) is the usual definition of 
creation and annihilation operators, the second condition 
is less standard; the motivation for this condition origi¬ 
nates from the fact that iLgff obeys 


[C„iLeff]=0, i = 


(One can verify this identity with straightforward alge¬ 
bra). As a result, we can restrict to a operators that 
commute with the CAs, and we will still have enough 
quantum numbers to completely diagonalize Hgg. 

We will call the a operators with E > 0, “annihilation 
operators”, and those with E < 0 “creation operators.” 
We will denote the annihilation operators by ai,...,aK, 
the creation operators by a\, ...,a^, and the correspond¬ 
ing E's by El,..., Ek {Ei > 0). 

The creation and annihilation operators have several 
important properties. To explain these properties, we 
need to recall some notation from section HA. There, 
we noted that there exists a change of variables C' = 
with V an integer matrix with determinant 
±1 and X defined by (15), such that the M x M matrix 
of commutators [C'',C'] takes the form 

/0i -V 0 \ 

[ClC'j]=2m\v 0i 0 , (B3) 

yo 0 Om-27/ 


with 


V 


/ di 

0 . 

. 0 \ 

0 

d2 . 

. 0 


0 . 

• djJ 


(B4) 


Here I is an integer with 0 < I < M/2 and 0/ denotes 
an I X I matrix of zeros, and similarly for Om- 2/- 





26 


With this notation, we can state the first property of 
the creation/annihilation operators: the number of lin¬ 
early independent annihilation operators is exactly 


K = N - M + I (B5) 

where I is defined as above. The second property of these 
operators is that they can always be chosen so that 

[ak,al,] = Skk' , [ak,ak'] = [a\,a\,] = 0 (B6) 

The third property of these operators (which 
is a consequence of the previous two) is that 
each of the {xi, operators 

can be written as a linear combination of 
{oi, aK^ o\,a\^, C2j^i, C'm^ ni,II m}. 

In appendix E we show that the above proper¬ 
ties are guaranteed to hold provided that we make 
one (rather technical) assumption: we assume that 
there is no operator R which is a linear combination 
of {xi, ...,Xn,Pi, ■■■,Pn}, is linearly independent from 
{Cl, ...,Cm}, and satisfies 


[R,H,s] = [R,C,] = 0 (B7) 

for all i. To understand the physical meaning of this as¬ 
sumption, note that if such an R operator existed, that 
would imply the existence of a continuous real-valued 
quantum number that we could use to label the low en¬ 
ergy eigenstates of H. Such continuous quantum num¬ 
bers cannot occur in finite-sized systems (and they don’t 
occur in any of the examples discussed in this paper) so 
we do not sacrifice much generality in making this as¬ 
sumption. 

Using these properties we can derive an important re¬ 
sult: the Hamiltonian Hgg can be written in the form 

K 

H,s = Y,Ekalak + F{C'2i+i,...,C'^) (B8) 

fc=i 

where F is some quadratic function. One way to prove 
this result is to observe that (i7eff — Sfc com¬ 

mutes with many other operators. For example, 

[ak' , ^ Eka\ak] = 0 (B9) 

k 

as one can see from the commutation relations between 
the ak operators and i/eff- Likewise, one can see that 

[a\,,Hea-'^Eka\ak\=Q (BIO) 

k 

At the same time, we have 

[C„ iLeff - ^ Ekalak] = 0 (Bll) 

k 

since = [Cj,ak\ = [Cj,al] = 0. To derive the 

consequences of these identities, we use the fact that the 


Xi,pi operators can be written as a linear combination of 
ak,a\,Ilj, and ■■■,C'j^. Clearly, this result implies 

that (iLeff — J2k ^koiak) can be written as a quadratic 
function of these operators. Examining the commutators 
(Bll) we see that this quadratic function cannot contain 
any H^- operators since these would make the commuta¬ 
tors with Cj nonzero. Similarly, from (B9) and (BIO), 
we can see that this quadratic function cannot contain 
any ak, a\ operators since these would make the com¬ 
mutators with and ak nonzero. We conclude that 

{Hes - J2k Ekalak) must depend only on ..., 

The expansion in Eq. B8 follows immediately. 


2. Step 2: Occupation number basis 


In the next step we construct a basis for the Hilbert 
space R. that is analogous to the conventional occupation 
number basis for a harmonic oscillator. To this end, we 
note that the following operators all commute with each 
other (see Eq. B3): 


(B12) 

Furthermore, these operators commute with the occupa¬ 
tion number operators {a|ai,..., }■ Therefore, we 

can simultaneously diagonalize (B12) along with {a^a/c}. 
We denote the simultaneous eigenstates by 


\9i, ...,0i,ipi, ...,ipi,x'j_^_i, ...,x'M-i,ni, ...,nK) 

or in more abbreviated form, \6, ip, x' , n) . Here the quan¬ 
tum numbers are defined by 

n) = e*®’ |0, ip, x', n), i = I, ..., / 

\0, p, x', n) = \9, p,x',n), i = l ^ 1,..., 2/ 

C[\e,p,x',n) = 2'Kx\_j\e,p,x',n), i = 2I+ 1,...,M 
alak\0,p,x' ,n) = nk\0,p,x',n), k = l,...,K (B13) 

where 0 < 9i,pi < 27r, while x'^ is real valued and 
ranges over non-negative integers. Importantly, one can 
show that there is exactly one simultaneous eigenstate 
\6, p, x', n) for each choice of 6, p, x', n so our labeling 
scheme is well-defined (see appendix F for a proof). By 
construction the \0,p,x',n) states form a complete or¬ 
thonormal basis for the Hilbert space H: these are the 
basis states that we seek. 

In fact, not only do the \9, p, x', n) states form a basis 
for R, but a subset of these states form a basis for the 
low energy Hilbert space ReS- To see this, recall that ReS 
consists of all states \ip) that satisfy cos{Ci)\tp) = \tp) for 
all i. Applying the two identities (A28-A29), we can see 
that ReS can be equivalently defined as the set of all 
states |^/>) satisfying 

cos(C'()|V’) = ... = cos(C'm)IV’) = \^) 


(B14) 


Substituting \9,n) into the above definition, we 
conclude that \6, ip, x', n) belongs to 'Heff if and only if 

9 = (27rai/di,27ra//d/), 

for some Ui = 0,1 ,di — 1, and 

¥> = (0, x' = {qi,...,qM-2i) 

with the qi’s being integers. It follows that the above 
states form a basis for Heg- We will denote these states 
using the abbreviated notation {|Q:,q,n)}, where a. and 
q are defined above, and n = (ni,..., uk) is the usual set 
of occupation numbers. 

3. Step 3: Eigenstates and energies 

We now have everything we need to diagonalize H^ff. 
Indeed, we’ve already seen that the \a,q,n) states form 
a basis for the low energy Hilbert space HeS- At the same 
time, from Eq. B8 we can see that the \cy.,q,n) states 
are eigenstates of HeS with energies 

K 

E{ol, q, n) = E UkEk + F{27Tqi,27TqM-2i) (B15) 

k^l 

We therefore have found all the eigenstates and energies 
of 

Now that we have the energy spectrum, we should 
point out one of its most important features: the en¬ 
ergy E is independent of the quantum numbers a. = 
(oi,..., a/). Since ai ranges from 0 < — I, it 

follows that every eigenvalue of ideff bas a degeneracy of 
at least 

I 

D = \{d, (B16) 

i=l 

Before concluding, there is one more issue we need to 
discuss: we haven’t yet explained how to determine the 
quadratic function E. One approach for computing F is 
to define a new set of operators $ 27 + 1 , •••, d’M by 

M K 

= '^WijIVj E'^{xikak +x*i,a\) (B17) 

3 = 1 k=l 

where the coefficients Wij and Xik are defined by Wij = 
and x^k = These coef¬ 

ficients have been chosen so that the operators have 
simple commutators with other key operators. In partic¬ 
ular, they have been chosen so that 

= 27ri<5,„ = [$„ 4] = 0 (B18) 

Thus, the ^ 21 + 1 , are conjugate variables to 

Once we have defined these operators, we can com¬ 
pute F by considering the commutators [<I>i,iJeff] and 


27 

ideff]]. Indeed, from the expression (B8), we de¬ 
rive 

dF 

= 2777 — 

^q^ 

d^F 

= (B19) 

These relations completely determine F since it is a 
quadratic function of the qi’s. (Alternatively, one can 
often compute F using problem-specific approaches, as 
in section III). 

4. Another method for finding creation and 
annihilation operators 

We have seen that one of the key steps in computing 
the energy spectrum of ideff is finding creation and anni¬ 
hilation operators a — that is, finding solutions to Eqs. 
BI,B2. Here we point out that there is an another, more 
convenient, way to formulate these equations. In this al¬ 
ternative formulation, we search for all operators a that 
obey 

M 

[ci, 77o] = Aa -f (B20) 

j=i 

[a,a]=0, i = (B21) 

for some scalars E, Xj with E 0. Here the Xj can be 
thought of as a kind of Lagrangian multiplier. We will 
show below that the above equations are mathematically 
equivalent to the previous equations (Bl) and (B2) in the 
sense that every solution to (B20) and (B21) is a solution 
to (BI) and (B2) and vice versa. The above equations 
are often more convenient than (Bl) and (B2) because 
they are written in terms of Flo and thus do not require 
us to compute Tfeff- 

To prove the equivalence between the two approaches, 
let us suppose that a obeys (B20) and (B21) for some 
E, Xj. We wish to show that a also obeys (Bl). To prove 
this, we take the commutator of (B20) with Q. The 
result is: 

M 

[a,[a,Ho]] = Y.Xj[a,[Cj,Ho]] 

7=1 

We deduce that 

1 “ 

1 ^ 

= (B22) 

i=l 

where in the second step, we used the Jacobi identity. 
Substituting the above formula for Xj into (B20), we de- 
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rive 

1 “ 

[a, Ho] + ^dCj,Ho] ■ [a, [Q,ifo]] = Ea 

i>*=i 

which can be rewritten as 

M 

[a, Ho] - ^ ■ [o,ni] = Ea 

j,i=i 

This is exactly (Bl). Conversely, one can check that any 
solution to (Bl), (B2) provides a solution to (B20) with 
Xj given by Eq. B22. 

Appendix C: Finite U corrections 

When U is large but finite, i7eff only gives the approxi¬ 
mate low energy spectrum of H. It is natural to wonder: 
what kinds of corrections do we need to add to i/eff if 
we want to obtain the exact effective theory at finite U1 
The goal of this appendix is to address this question. 

1. Simple example 

We begin by discussing the example from section II A: 

H = ^ -\- ^ - U cos(27ra;). For this example, the 

low energy effective Hamiltonian in the infinite U limit 
is Hes = —f- while the low energy Hilbert space HeS 
is spanned by position eigenstates Hg)} where q is an 
integer. 

To understand the finite U case, let us imagine re¬ 
peating the derivation of Heff from appendix Ala, but 
without taking the limit C/ —>■ oo or making any other ap¬ 
proximations. In such a hypothetical exact calculation, 
we would first write H as a sum H = Hi -\- H 2 where 
Hi = — U cos(27ra;) and H 2 = Next, we would 

find find the exact eigenstates and energies of Hi. Ac¬ 
cording to Bloch’s theorem, these states can be labeled 
as \k, n) where n = 0,1,... is the band index and k is the 
crystal momentum, — tt < k < tt. To complete the calcu¬ 
lation, we would construct a low energy effective theory 
describing the energy spectrum of H below the first band 
gap by treating H 2 as a perturbation to Hi , and including 
terms to all orders in perturbation theory. In this way, we 
can imagine deriving an effective Hamiltonian H^ff giving 
the exact low energy spectrum of H in the finite U case. 
This effective Hamiltonian i?eff would be defined in an 
effective Hilbert space HeS spanned by the states in the 
lowest band, {|A:,0)}. Equivalently, we can describe the 
effective Hilbert space, T-LeS as the span of the Wannier 
states {Ig)} defined by |g) = 0)dfc. Here, Ig) 

denotes the Wannier state localized near the cosine min¬ 
imum at a; = g. 

Although we will not perform the above calculation, 
we can still make some qualitative statements about the 


structure of the final result. Indeed, from our analysis of 
the infinite U limit, we know that the resulting i7eff must 
take the form 

H.S = E ^l9)('?l + (CI) 

q qq' 

where —>■ 0 as ?7 —>■ 00 . Furthermore, we can esti¬ 
mate the size of e^g' for finite U. There are two cases to 
consider: q ^ q' and q = q'. When g 7 ^ g', we expect that 
^ 99 ' const.-^/mU gjjj^Qg these off-diagonal terms are 

generated by tunneling between different cosine minima 
of Hi. As for the q = q' case, these terms can be esti- 
mated as Cgg ~ jj since the leading order contribution 
to Cgg comes from the second order terms in the pertur¬ 
bative expansion in H 2 ^ computed in Eq. A4. 

It is instructive to rewrite the above expression (CI) 
for Hes in terms of the operators x,p. In this language, 
we have 

H,s= — + E (C2) 

n— — oo 

where e„ is a function satisfying en(q') = f-{q'-n)q' for all 
integers g'. Here, the equivalence between (CI) and (C2) 
follows from the fact that e™P|g') = |g' — n). Translating 
our results about Cgg' into this alternative language, we 
see that e„ has an exponential dependence on t/U for 
n 7 ^ 0, and scales like 1/U for n = 0. 

Equations (CI, C2) give the qualitative structure of 
i?eff in the finite U case. To obtain more quantitative 
results, we would need to explicitly compute the coeffi¬ 
cients Cgg' or the function en{x) as a function of U, K, m. 
In principle it should be possible to compute these quan¬ 
tities in a large U expansion — for example using the 
instanton approach outlined in section 7.2.3 of Ref. 35. 

2. General case 

We now discuss the finite U corrections for more gen¬ 
eral systems. First, we consider a scenario in which 
only one of the cosine terms has a finite coefficient while 
the others are taken to be infinitely large. In other 
words, we consider Hamiltonians of the form H — Ho — 
U cos(C'i) — U' cos(Cj) in the limit where U is finite 
but U' —>■ 00 . In this case, the finite C/ corrections only 
generate tunneling processes of the form Ci —>■ Ci — 2?™; 
other tunneling processes, Cj —>■ Cj — 27rn, are suppressed 
by t/'. It follows that the finite U corrections must com¬ 
mute with {Cl, ...,Ci-i,Ci+i, ...,Cm} but don’t have to 
commute with Ci. At the same time we know that the 
most general low energy operator is of the form shown in 
Eq. DIO. Combining these two facts, we conclude that 
the finite U corrections can be written in the form 

00 

E e“n-e„({afe,4,C',+J) (C3) 
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for some functions e„(ai,a^, a|,a^, C 2 /+ 1 ',C^)- 
Note that these terms are generalizations of the finite U 
corrections (C2) that we discussed for the above example. 
As in the example, the functions have a U dependence 
which we have not shown explicitly. In particular, the e„ 
with n ^ 0 vanish exponentially in '/U as C/ —>■ 00 since 
these terms are generated by (non-perturbative) instan- 
ton effects, while cq vanishes like 1/C/ since this term 
originates from perturbative corrections. 

Next we consider the case where all the cosine terms 
have finite coefficients. In this case, all tunneling pro¬ 
cesses Ci ^ Ci — 2'Kmi are allowed so the finite U cor¬ 
rections take the form 

. e^({afc,4,C',+J) (C4) 

m 

with the sum running over M component integer vectors 
m = {mi, As above, the e™, are unknown func¬ 

tions of {oi,..., Ofc, 4 ,..., 4 , (727+1',..., which also 
depend on U. The Cm with m ^ 0 vanish exponentially 
in \/U as C/ —>■ 00 , while eo vanishes like 1/C/. In princi¬ 
ple, it should be possible to compute Cm using instanton 
methods^®, but we won’t discuss this computation here. 


3. Splitting of ground state degeneracy 

One application of this formalism is that we can use 
it to analyze how the H-fold ground state degeneracy of 
i/eff splits at finite U. Indeed, according to lowest or¬ 
der perturbation theory, we can determine the splitting 
of the ground state degeneracy by projecting the finite 
U corrections onto the ground state subspace and then 
diagonalizing the resulting D x D matrix. This diago- 
nalization problem is system dependent so we cannot say 
much about it in general, but we would like to mention a 
result that is useful for setting up the computation. This 
result applies to any system in which the commutator 
matrix Zij is non-degenerate, i.e. it applies to the case 
where M = 21. The result states that the matrix ele¬ 
ments of the finite U corrections are proportional to the 
matrix elements of (where Tj is defined as in 

Eq. 24): that is, 

where |Q;),|a') are ground states and Um is some un¬ 
known proportionality constant. Although this relation 
does not tell us the value of Um, it is still useful since it 
tells us the form of the matrix that we need to diagonal¬ 
ize. 

To derive equation (C5), consider the difference 11^ — 
Tj . This difference is a linear function of x^’s and p^’s 
so we know it can be written as a linear combination of 
{afc,aj,,(7'} as in Eq. D5.^° At the same time, it is easy 
to see that II^- — T^ commutes with the Ci operators and 
hence also the C{ operators. It follows that the expansion 


of Hj — Tj in terms of {ofc, a\, cannot contain any C[ 
operators: that is, 

Ilj = T^- -I- (linear combination of Cfc and a\). 

If we exponentiate this relation, Eq. C5 follows easily. 
Note that the constant prefactor Um comes from the 
and a\ operators appearing in this expression. 


Appendix D: Matrix elements and operators 
1. Matrix elements 

If we wish to have a complete low energy theory, we 
need to do more than just find the energies of the low- 
lying states: we also need to be able to compute matrix 
elements of operators between these states. To address 
this issue, we now outline a general procedure for com¬ 
puting matrix elements (a, q, n\0\cx, q, n). 

Our basic strategy is as follows: suppose we are 
given a general operator O which is some function of 
{xi, ...,Xn,Pi, ■■■,Pn}- What we will do is express O as 
a function of the operators ak,a\,Ci,Ili, etc. Then we 
will use the known matrix elements of a^, a|,, (7+ I!,, etc 
to compute the matrix elements of O. 

We now demonstrate how this approach works in more 
detail. We focus on two cases which are particularly 
important: {!) O = A and (2) O — where A is 
some linear combination of Xi’s and pi’s. We begin 
with the first case, O = A. In this case, our strat¬ 
egy will be to express 7l as a linear combination of 
{oi,..., otc, a\,C 2 j.^_i, ■■■, C’mj Hi, .., IIm}: 

K MM 

= ^^(KfeU/c + AfeoJ.) + PiCl (Dl) 

fc=l i=27+l i=l 

We know such an expression must exist, since the x^’s 
and Pi’s can be written as linear combinations of these 
operators, as shown in appendix B 1. 

Our task is now to find the expansion coefficients, 
Kfc) Afc, Pi, t'i. This can be accomplished by consider¬ 
ing appropriate commutators. For example, if we take 
the commutator of Eq. Dl with Cj, we derive Vj = 
— ^ [Cj , A] . Similarly, taking the commutator with a\ 
and Ofc, we derive 

M 

Kk = -[al,A]+^v,[al,T\,] 

M 

[^k 1 ^ ^ \^k •) n^] 
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Finally, taking the commutator with Ilj, we hnd 

M K 

— [nil'll ^ ^ 

K M 

— Afc [Ilj, a^]) — [Ilj, Hi] 

k—l i—1 

which we can invert to find ^i. 

Once we have the expansion (Dl), our problem reduces 
to finding the matrix elements of the operators Ofc, aj,, 11 ^ 
and C' 2 j_|_]^,..., C^. The matrix elements for a^, 0 ^,( 7 ' 
can be written down without any work: 

afc|a,q,n) = y^ja, q, n - e*) 

a\\oL,q,n) = y/nu + l|a,q,n + ei) (D2) 

C'-|a,q,n) = 27rqi_2/|Q:, q, n), i = 2/ + 1, ...,M(D3) 

Here e* denotes the if-component vector ej = 
(0,..., 1 ,0) with a “ 1 ” in the ith entry and 0 every¬ 
where else; thus the first two equations encode the fact 
that the Ofe, aj, act as raising and lower operators in the 
occupation numbers n. As for the last equation, this 
follows from the definition of |a, q, n) in section IIB. 

All that remains are the H^ operators. We now ar¬ 
gue that the matrix elements of these operators vanish 
identically:^^ 

{a.',q',n'\I{i\a.,q,n) =0 (D4) 

The derivation of Eq. D4 follows from two observa¬ 
tions. First, we recall from our derivation of iJeff that 
the low energy states \a.,q,n) are degenerate ground 

states of the Hamiltonian Hi = ~ 

UJ2i=i‘^o^iCi) in the limit U ^ 00 . Thus, we know 
that Hi ja, q, n) = E\a., q, n) for some scalar E. Second, 
we note that H^ = ^ Hi], so in particu¬ 

lar Hi is a linear combination of commutators [Cj,Hi]. 
Equation (D4) now follows from the identity 

{a',q',n'\[Cj,Hi]\ct,q,n) = 

(a', q', n'\CjE — ECj\a, q,n) = 0 

Putting together equations (D1,D2,D3,D4), our task 
is complete: we can compute matrix elements of 
any operator A that is a linear combination of 
, ..., Xtv : Pi, •", Pn'\ • 

We now move on to the second case: opera¬ 
tors of the form where is a linear combi¬ 

nation of Xi’s and pi’s. For this case, our strat¬ 
egy is to express A as a linear combination of the 
operators {oi,..., a^f, a|,..., ,..., C^} along with 

$ 2 /-i-i,..., (where the operators are defined in Eq. 
B17): 

K MM 

A = '^{Kkak + Xkol) + y] Pi^i + y] ViCl (D5) 

k—l 4=27+1 2=1 


We know such an expression must exist because (1) the 
above set contains 2N operators, as can be seen from Eq. 
B5, and (2) these operators are linearly independent, as 
one can show using the same arguments as in appendix 

E. 

Just as in the previous case, we can find the ex¬ 
pansion coefficients Kk, Xk, Pi,k'i by taking appropri¬ 
ate commutators. After finding the coefficients, our 
problem reduces to computing the matrix elements 
of The matrix elements of 

^ found straightforwardly using (D2). 

As for one can show, using the commutation rela¬ 

tions (B18), that acts as 

e^"'^'\a,q,n) = \a,q±ei_ 2 i,n) (D 6 ) 

for i = 21 + 1,...,M. Thus, act as raising and 

lowering operators on the q quantum numbers. 

Finally, we need to discuss the matrix elements of 
There are three cases, each of which needs to 
be treated differently: ( 1 ) 1 < * < /, ( 2 ) /-|-1 < i <21, 
and (3) 2/ -I- 1 < i < M. For the first case, Eq. B13 
implies that 

+, q, n) = |q,^ z = l,...,/ (D7) 

For the second case, one can show using the commutation 
relations [C'',C''] = 27 riZb that 

|q,^ n) = \oL± Si^i, q,n), i = I + 1 ,..., 21 

(D 8 ) 

where the addition is performed modulo di. Note that 
equations (D7-D8) imply that the operators act 

like “clock” matrices for i = 1 ,...,/, while the operators 
^±iCJdi_j “shift” matrices for i = I + 1,...,2/; 

thus these operators generate a generalized Pauli algebra. 
Finally, for the third case, the matrix elements of 
can be obtained by exponentiating Eq. D3. 

The above equations tell us everything we need to com¬ 
pute (a', q', n'|e*-^|Q;, q, n) for the special case where the 
coefficient pi in (D5) is an integer for every i, and Vi is 
an integer multiple of 1/di for i = 1,...,I, and a multi¬ 
ple of 1/di-i for i = / -I- 1,..., 2/. To complete the story 
we need to explain how to evaluate matrix elements if 
the Pi and Ui coefficients aren’t quantized in this man¬ 
ner. To understand this case, suppose that pi is not an 
integer. It then follows that the commutator [A, C/\ is 
not an integer multiple of 27ri, which implies that the 
state q, n) is an eigenstate of C' with non-integer 

eigenvalue. But this means that e*^|Q:,q,n) is orthogo¬ 
nal to the low energy Hilbert space, Hee so that all the 
matrix elements (a', q', q, n) vanish identically. 

Similar reasoning shows that these matrix elements also 
vanish identically if the Vi coefficients aren’t quantized as 
above. Hence, the above special case is actually the only 
case where the matrix elements are nonzero. 
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2. Low energy operators 


We now turn to the question of how to describe the 
most general operators in the low energy theory — that 
is, the most general operators that act within the low 
energy Hilbert space ’HeS- The simplest way to do this 
is to write the operators as linear combinations of outer 
products of the form \ol', q',n'){a.,q,n\. This approach 
is straightforward, but it can be unwieldy since many 
operators look complicated in this representation. 

Alternatively, we can represent any operator O in the 
low energy theory as an infinite power series in the oper¬ 
ators 


{C'27+ii ■••I 6 


Ki 

±i^ 2 /+i 








That is. 




(D9) 

for some function / that can be expanded as a power 
series. Indeed, to prove this result, we need to establish 
two facts: ( 1 ) we need to show that the above operator 
O maps ’Heff —'He®, and ( 2 ) we need to show that the 
functional form for O is sufficiently general that it can 
reproduce any linear map from "Heff 'HeS- The first 
fact is easy to prove: we can see that all of the operators 
{afc, 4 ,C''j+,, commute 

with {cos(Ci),..., cos(Cm)} and therefore O also has this 
property. The second fact is harder to prove, but can be 
established straightforwardly by examining the matrix 
elements for 

given in equations (D2, D3,D6,D7,D8). 

To get more intuition about the above 
representation, we note that the operators 


have a nat¬ 
ural interpretation in terms of the phase space of the 


low energy theory: the {afc,a 4 operators are conjugate 
variables that describe real degrees of freedom at low 
energies; the operators describe angular-valued 

degrees of freedom and the {C' 2 /_|_j} describe the corre¬ 
sponding conjugate discrete degrees of freedom; finally. 


the operators are conjugate variables 

that describe the finite discrete degrees of freedom. The 
latter operators can be thought of as generalized Pauli 
operators, similar to 

In fact, there is yet another way to parameterize the 
low energy operators which is often more convenient to 
use: every low energy operator O can be written as 
an infinite power series in {ak,al.}, {^ 27 + 1 ,..., and 
That is. 


0 = /({afe,4>^W»e±®-}) 


(DIO) 


To derive this parameterization, it sufhces to show that 
we can express the operators and 


{g±jCj^i/di} fashion; once we establish this prop¬ 

erty then this parameterization follows immediately from 
the previous one (D9). To derive the latter property, we 
use the definition of (BI7) which states that <I>i is a 
linear combination of the Ilj,ak,al. operators. Impor¬ 
tantly, if we examine this expansion, we can see that the 
coefficients of the H^’s are integers. Therefore if we ex¬ 
ponentiate (BI7), we immediately see that can be 
written as a monomial in multiplying 

a power series in {a^, a^- Identical reasoning shows that 
^±iCJdi g±iCjg.;/<ii ggg gjgQ j^g written in this way; 
this completes the proof. 


3. Ground state operators 

So far we have focused on operators acting within the 
low energy Hilbert space HeS- However, in some cases we 
may be interested in even lower energy scales in which 
we will need to think about the D-dimensional subspace 
spanned by the set of degenerate ground states. In our 
notation, these ground states take the form | q ;, q, n) with 
n — q = 0. We will label them by |a) = |a, 0, 0). 

Once we consider the ground state subspace, we face a 
similar question: how do we describe the operators that 
act within this subspace? The simplest way to do this 
is to write the operators as linear combinations of outer 
products, | q :')( q :| but this representation is not always 
the most convenient one. 

Another way to represent the operators O that act 
within the ground state subspace is to write them as 
polynomials in 

|g±»c;/di ^ g±jc;/dj ^ g±*c;+i/<ii ^ g±*cb/d7 

That is, 

O =(DII) 

for some polynomial function /. Indeed, to prove this 
result, we need to establish two facts: (I) we need 
to show that the above operator O maps the ground 
state subspace to itself, and ( 2 ) we need to show that 
the functional form for O is sufficiently general that it 
can reproduce any linear map from the ground state 
subspace to itself. To prove the first fact we note 
that the operators commute with 

{cos(C'i),..., cos(Cm), Heff} and therefore O also has this 
property. The second fact is harder to prove but can be 
established straightforwardly using the matrix elements 
for given in Eqs. D7,D8. 

There is also a third way to parameterize the ground 
state subspace operators which applies to the case 
where the commutator matrix Zij = 2 iy[Ci)Cj] is non¬ 
degenerate. This alternative parameterization involves 
the operators Ti, ...,rM, defined in Eq. 24. More specif¬ 
ically, this parameterization states that general ground 
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state operators O can be written as polynomials in 
That is, 

0 = /({e±*r.}) 

for some polynomial /. To derive this result, it suffices 
to show that we can write and as such 

polynomials; once we establish this fact, the parameteri¬ 
zation follows immediately from the previous one (Dll). 
To prove the latter property, we note that the opera¬ 
tors and can be equivalently written 

as exp(±f We then use the following 

identity, which can be derived easily from the definition 

ofq,r,: 

M M 

1=1 1=1 

The important point about this identity is that {V~^)ji 
is an integer matrix, so we deduce that the left hand side 
can be written as a linear combination of Tj’s with integer 
coefficients. It follows that exp(±i can 

be written as a monomial in which is of 

course a special kind of polynomial. This completes the 
proof. 

Appendix E: Properties of creation and annihilation 
operators 

In this appendix, we derive three properties of cre¬ 
ation/annihilation operators that were used in appendix 
B 1. These properties are guaranteed hold as long as we 
assume that there is no operator R which is a linear com¬ 
bination of {a;i, ...,xn,pi, ...,pn}, is linearly independent 
from {Cl, ...,Cm}, and satisfies 

[R, Cl] = ... = [R, Cm] = [R, H,s] = 0 (El) 

(See appendix B 1 for more discussion about this assump¬ 
tion). 

The first property that we will derive is that the num¬ 
ber K of linearly independent annihilation operators is 
exactly 

K = N - M + I (E2) 

We prove this result by establishing two opposing in¬ 
equalities: K < N — M + I and K > N — M + I. We 
start by showing the inequality K < N — M + I. 

To begin, consider the set of all operators that are lin¬ 
ear combinations oi {xi, ...,Xn,Pi, ■■■,Pn} and that com¬ 
mute with {Cl,...,C m}- This set forms a vector space 
since it is closed under addition and scalar multiplication. 
We will call this vector space V. 

Next, observe that if A is an operator that belongs 
to V, then the commutator [A,Heg] also belongs to V: 
to see this, it suffices to show that [Ci , A] = 0 implies 


[Ci, [A, iJeffj] = 0. The latter result is a simple applica¬ 
tion of the Jacobi identity: 

[C„ [A, iJeffj] = -[A, [J/eff, C,]] - [iJeff, [Q, A]] = 0 

where the second equality follows from the fact that 
[iJeff, Ci] = 0. In view of this property, we can think 
of the commutator map A —>■ [A, i/eff] as defining a lin¬ 
ear mapping from V ^ V. We will denote this linear 
mapping by S. 

In this language, the annihilation operators a corre¬ 
spond to eigenvectors of S with positive eigenvalue. Our 
task is thus to bound the number of these eigenvectors. 
To this end, we note that the total dimension of V is 
2N — M. At the same time, we know that the operators 
C 2 /_|_i,..., Cm belong to the space V and are eigenvectors 
of S with eigenvalue 0. It follows by dimension counting 
that S has at most (21V — M) — (M — 2I) = 2N — 2M + 21 
linearly independent eigenvectors with nonzero eigenval¬ 
ues. 

To complete the derivation of the first inequality, we 
note that the eigenvectors of S with nonzero eigenvalues 
come in ±C pairs: if [A, i/eff] = E ■ A then [Al, i/eff] = 
—//• Ab It follows that S has at most N — M + I linearly 
independent eigenvectors with positive eigenvalues. In 
other words, K < N — M + I. 

Now we prove the second inequality, i.e. K > N — M + 
I. The first step is to observe that physical considerations 
guarantee that i/eff is positive semi-definite. It follows 
that we can write i/eff as 

i^eff = E^’ (E3) 

i=l 

where each Bi is a real linear combination of 
{xi, ...,xm,Pi, ■■■,Pn}^ and where each Bi is linearly in¬ 
dependent from the others. Here L is some integer with 
L<N. 

Next, we observe that the Bi operators have the prop¬ 
erty that [Bi,Cj] = 0 for all i,j: one way to de¬ 
rive this fact is to expand the commutator [Cj , i/eff] as 
^^i[Cj,i/i]i/i, and then use the fact that the Bi op¬ 
erators are linearly independent as well as the fact that 
[C„//eff]=0. 

To proceed further, consider the L x L matrix = 
[Bj,Bi]. This matrix is imaginary and skew-symmetric, 
which implies that it is diagonalizable. Thus, we can 
find L linearly independent eigenvectors vi,...,vl with 
eigenvalues Ei,E^. The corresponding operators = 
obey [fli, i/eff] = Eitti- Furthermore, we know 
that [ai,Cj\ = 0 since [Bi,Cj] = 0 for all i,j. 

Clearly the Oi operators obey almost all of the con¬ 
ditions for annihilation operators. All we have to show 
is that we can find a subset oi N — M + I linearly in¬ 
dependent tti’s with eigenvalues Ei > 0. To this end, 
consider the set of all operators that are linear combi¬ 
nations of {xi, ...,xn,Pi, ...,piy}, and commute with all 
the Ci and operators. This set is a vector space since 
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it is closed under addition and scalar multiplication. We 
will call this vector space W. Let us try to find the di¬ 
mension of W. To do this, note that any operator in W 
commutes with all the Bi operators, since the Bi oper¬ 
ators can be written as linear combinations of the Oi’s. 
It then follows that any operator in W must also com¬ 
mute with ffeff- But by our assumption (El), the only 
operators that commute with both ffgg and the Ci op¬ 
erators can be written as linear combinations of Ci, or 
equivalently, linear combinations of C)^^}. We 

conclude that the dimension of W is at most M — 21. 

Given that the dimension of W is at most M — 21, it 
follows that there must be at least 2N — M + 21 linearly 
independent Oi, Ci operators. Hence, at least 2N — 2M -|- 
21 of the tti operators are linearly independent from the 
Ci operators. It then follows from the assumption (El) 
that at least 2N — 2M + 21 of the Oi’s have Ei ^ 0. At 
the same time, we know that the Ei eigenvalues come 
in pairs of opposite sign since Wj is skew-symmetric and 
imaginary, so we conclude that at least N — M + I of 
the Oi’s have positive eigenvalue. This establishes that 
K > N — M + I and completes the proof of Eq. E2. 

We now prove the second property of the creation and 
annihilation operators. This property states that the aCs 
can always be chosen so that 

[ak,a\,] = 5kk', [ak,ak'] = = 0 (E4) 

We begin with the first relation. To prove this, we ob¬ 
serve that [ak,a\,] = 0 unless Ek = E(,. Indeed, this 
follows from the Jacobi identity: 

\(lk , \^k' ’ 7 Fife]] [Ueff? \^k ; 

Now hx a particular Ek- There are two cases to con¬ 
sider: Ek may be degenerate or non-degenerate. If Ek is 
non-degenerate then we can simply normalize the corre¬ 
sponding ttfe such that [afc,aj)] = 1; the vanishing of the 
commutators between ak and other a^,’s is guaranteed. 
On the other hand if Ek has some degeneracy, then we 
can use the Gram-Schmidt procedure to choose the asso¬ 
ciated Ofe’s so that they obey [a/c,a^,] = Skk'- Again, all 
the other commutators vanish automatically. Thus, in all 
cases we can choose the operators to satisfy the first 
relation. As for the second relation, this follows from the 
Jacobi identity by similar reasoning. 

Einally, we move on to prove the third property 
of the creation and annihilation operators. This 
property states that each of the {xi, ...,xn,pi, ...,pn} 
operators can be written as a linear combination 
of (ui,..., ,..., , G2,..., , H]^ ,. Hjvf }. To 

prove this, we only need to show that the operators in 
this set are linearly independent since (I) there are 2N 
of them all together according to Eq. E2, and (2) they 
are all linear combinations of {xi, ...,X]sf,pi, ...,pn}. The 
fact that they are linearly independent can be seen as fol¬ 


lows: suppose that 

K MM 

'y ( (/tfeOfe -i- Afco^) -T y ( +iCi + y ( i^iHi = o 

fc=l i=2I+l i=l 

for some scalars Kk, Xk, Then, if we take the com¬ 

mutator of both sides with Cj , we immediately see that 
i^j = 0 since [ak,Cj] = [a\,Cj] = [C[,Cj] = 0 while 
[Hi, Cj] = Sij. We therefore have 

K M 

'^{Kkttk + Xkal) + ^ piCl = 0 

fe=l i=2I+l 

Next, we take the commutator of both sides with and 
we deduce that = 0 since [ak, G'] = 0 while [uk, a^,] = 
Skk'- Similarly, taking the commutator with aj, shows 
that Kk = 0. Our relation now becomes 

M 

E = 0 

1=27+1 

Einally we note that the G' are all linearly independent 
by construction, so all the pi^s must vanish. Hence, all 
the coefficients in our original expansion must vanish, 
which implies that the above operators are linearly inde¬ 
pendent, as claimed. 

Appendix F: Uniqueness of simultaneous eigenstates 

In this appendix, we show that there is exactly one 
simultaneous eigenstate \6,ip,x',n) for each choice of 
6,cp,x',n. Here 9,4> are I component angular-valued 
vectors, while n is a AT component non-negative integer 
vector, and x' is a M — 21 component real valued vector. 

To see that there is at least one eigenstate for each 
choice of these quantum numbers, we observe that the 
explicit formula in Eq. A30 implies the weaker result 
that there is at least one state for each choice of 9, (p, x'. 
Now consider the subspace of states with hxed values 
oi 9, (p,x'. It is clear that the creation and annihilation 
operators a^, ak map this subspace onto itself. Therefore, 
using the same arguments as in the algebraic analysis 
of the harmonic oscillators, it is easy to see that this 
subspace contains at least one state for each choice of 
occupation number n. 

Gonversely, to see that there is at most one state 
\9, (f, x', n) for each choice of 9, p, x', n, we recall that all 
the Xi,pi operators can be expressed as a linear combina¬ 
tion of {ai,...,aif, 4 ,...,a^,G 2 j+i,...,GV,ni,..,nM}- 
One application of these expressions is that we can use 
them to compute the expectation value of any opera¬ 
tor of the form = exp(z -I- XiPi)) in 

any state with quantum numbers 9,ip,x',n. This com¬ 
putation is completely algebraic and depends only on 
the quantum numbers 9,ip,x',n as well as the parame¬ 
ters K, A. Thus we conclude that the quantum numbers 
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9,(p,x',n completely fix the expectation values of the 
Ok„x operators. But the above operators Ok,,a are suf¬ 
ficiently general that any operator f{xi,Pi) can be con¬ 
structed by taking an appropriate linear combination of 
them. Hence, if two states share the same quantum num¬ 
bers 0,cp,x',n, then they must have the same expecta¬ 
tion values with respect to all operators in the Hilbert 
space and hence must be equivalent to one another up to 
a phase. 


Appendix G: Regularizing the cosine term 

In this appendix, we revisit the problem of a fractional 
quantum spin Hall edge with a single magnetic impurity: 

H = Ho-Ucos{C) (Gl) 

The new element in our discussion is that we regularize 
the argument of the cosine term, replacing C = k((j)^{0) + 
</'t(0)) with 

/•V2 

C= k{(j)^{x) + (pi{x))6{x)dx 

J-L/2 

where <5 is an approximation to a delta function — i.e. a 
narrowly peaked function with f 5{x)dx = 1. After mak¬ 
ing this replacement, we repeat the analysis in section 
HIB in which we constructed creation and annihilation 
operators for the low energy effective Hamiltonian Ideff- 
Our main result is that we find that when the cosine term 
is properly regularized, the condition [a, C] = 0 translates 
to the constraint in Eq. 41; otherwise the regularization 
doesn’t change much. 

As in section HIB, our task is to find all a operators 
such that (1) a is a linear combination of our fundamental 
phase space operators ^^(yo), </>t(?/o)} and 

(2) a obeys 

[a, Hq\ = Ea + A[C', Hq\ + Hq] + A^[(5^) Hq] 

[a,C] = [a,Qt] = [a,Qt]=0 (G2) 

for some scalars E, A, A^, A^ with E ^ 0. 

Given that [a, Q^] = [a, = 0, we deduce that 

4’^(yo)i4‘i(.yo) cannot appear in the expression for a. 
Hence, a can be written in the general form 


get 


f^{y) = - Q^{y)) + A^Q^iy)] (G4) 

h{y) = - 0;(y)) + i?20;(2/)] (G5) 

where p = E/v and 

A2 = Ai- Afc, B2 = Bi- Xk (G6) 

and 

0 ^(y) = r e-^P^5{x)dx 
J-L/2 

0 t(2/) = r e^^-S{x)dx 

J-L/2 

Note that both 0^(a;), 0^(a;) reduce to the Heaviside step 
function 0(a;) in the limit that S{x) is a delta function. 
Eliminating A from (G6) we see that 

A 2 — Ai = B 2 — Bi (G'^) 


So far the regularization hasn’t taught us anything 
new: all of our results are similar to what we found in 
equations (36,37,38), when we analyzed the unregular¬ 
ized cosine term with C = k{(j)^{0) + ^i^(O)). So why is 
the regularization important? The reason is that it al¬ 
lows us to properly treat the constraint [a, C] = 0, as 
we now demonstrate. The first step is to use (G3) to 
translate the constraint to 



Iftiy) - fdy)]hy)dy = o 


Next, substituting our expressions for /^, (G5), we de¬ 

rive 



- B,e-^Py)~Siy)dy 

- A,)EPye^{y))Siy)dy 
Bi)e-^pye^{y)yS{y)dy = 0 


(G8) 


/■i/2 

a = / [f/iy)c>vyiy) + fi{y)dy(l)yy)]dy (G3) 
J-L/2 

Next, from the first line of Eq. G2 we obtain the following 
set of differential equations, 

-iyf/{y) = Ef^iy) + \kiv5{y) 

= Efiiy) - Xkiv6{y) 

(The A^,A^ terms drop out of these equations since 
commute with Hq). Solving these equations, we 


While the constraint (G8) looks complicated, it simplifies 
considerably in the low energy, long wavelength limit, 
i.e. the limit where pb 1 where b is the width of 
the function S. In this limit, the first integral in (G8) 
evaluates to {Ai — Hi), but the second integral is a bit 
trickier. To compute this integral, we use the identity 

hm / 5{y)Qyy)e^^Py = - (G9) 

J_I^I 2 2 

where a =t, 4-. This identity can proved for the case a =t 
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by noting that 


the form 


fLj’l 


lim / 

pb^0j_[^/2 

ry 

= lim / / 

pb^0j_j^/2 J-L/2 


~6{y)Q^{y)e^^Pydy 

S{y)Six)e-^P^^^Pydxdy 


fL/^ rv ^ ^ 

= lim / / 5{y)5{x)dxdy 

pb^0j_i^/2 J-L/2 

1 ri/2 ri/2 

= lim - / / S{y)d{x)dxdy 

pfc-i-o ^ J-L/2 J-L/2 


1 

2 


The proof for the case cr =4, is similar. 

Applying the above identity (G9), the second integral 
in (G 8 ) evaluates to ^[(A 2 — Ax) — (i ?2 — Bi)] so that 
(G 8 ) becomes 


(^1 


, (Aa - Ai) - (B 2 - Bx) 

Bx) + ^ = 0 


N 

H = Ho- U'^cos{C,) 
2 = 1 


Ci = 2k{(j)^{xi) + (fb/ixi)) 


(HI) 

(H2) 


where Hq is defined in Eq. 31. Notice the factor of 2 in 
the definition of Ci. 

We can see that this system is identical to the mag¬ 
netic impurity system studied in section III C except for 
the above factor of 2 . As a result, almost all of the results 
derived in section IIIG carry over to this case without 
change. In fact, the only difference comes when we com¬ 
pute the commutator matrix Zij = 2 ^[Gi,Cj]. For the 
above case, takes the form 


0 2 - 


Zij — 




0 ••• 0 2-2 
-2 ••• -2 0 0 
\ 2 ••• 2 0 0 / 

In order to transform the commutator matrix into 
canonical form, we make the change of variables 


Simplifying, we obtain 


C7( = Cl, C 2 = C 3 = 27rQ^ -b 27rQ; 

C'^ = Cm-2 — Cm-3j 771 = 4, ..., N J- 2 


Ai + A 2 _ Bx + B 2 

2 “ 2 


This is precisely the constraint from Eq. 41. 


(GIO) 


In this basis we obtain 


/O -2 0 

'200 
0 0 0 




Appendix H: Degeneracy from spontaneously 
broken time-reversal symmetry 

In this appendix, we consider a FQSH edge in a disk ge¬ 
ometry with N time-reversal invariant impurities located 
at positions Xx, ...jXn. We show that when U —)■ 00 , the 
ground state is two-fold degenerate, which is consistent 
with spontaneous time-reversal symmetry breaking. 

Because we consider time reversal-invariant impurities, 
the dominant scattering process in this system involves 
two-particle backscattering. Thus, the Hamiltonian takes 


Vo 0 0 ••• 0/ 

We can see that ■ is in the canonical skew-normal form 
shown in Eq. 14, with the parameters M = N-\-2, I = 
dx = 2 . 

If we now compute the degeneracy using Eq. 20, we 
obtain D — dx = 2. We conclude that the ground state 
is two-fold degenerate (in fact, every state shares this 
degeneracy). This degeneracy makes perfect sense phys¬ 
ically since we expect that the edge will exhibit spon¬ 
taneous time-reversal symmetry breaking in the limit of 
large and such symmetry breaking naturally 

leads to a two-fold degeneracy. 
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